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We demonstrate that partonic correlations that would traditionally be identified as subleading 
on the basis of a generalized TMD-factorization conjecture can become leading-power because of 
TMD-factorization breaking that arises in hadron-hadron collisions with large transverse momentum 
back-to-back hadrons produced in the final state. General forms of TMD-factorization fail for 
such processes because of a previously noted incompatibility between the requirements for TMD- 
factorization and the Ward identities of non-Abelian gauge theories. We first review the basic steps 
for factorizing the gluon distribution and then show that a conflict between TMD-factorization and 
the non-Abelian Ward identity arises already at the level of a single extra soft or collinear gluon when 
the partonic subprocess involves a TMD gluon distribution. Next we show that the resulting TMD- 

t-H factorization violating effects produce leading-power final state spin asymmetries that would be 

classified as subleading in a generalized TMD-factorization framework. We argue that similar extra 

TMD-factorization breaking effects may be necessary to explain a range of open phenomenological 

J_( QCD puzzles. The potential to observe extra transverse spin or azimuthal asymmetries in future 

experiments is highlighted as their discovery may indicate an influence from novel and unexpected 

•^£ large distance parton correlations. 
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Qh I. INTRODUCTION 

&* . 

This paper examines the consequences of factorization breaking in inclusive high energy cross sections that are 

differential in the transverse momentum of produced particles, with a focus on the region of small transverse mo- 
mentum where intrinsic motion associated with hadron structure becomes significant. There is at present a wide 
i-H range of motivations for predicting and measuring intrinsic transverse momentum effects experimentally, so much 

^ effort continues to be devoted to developing methods to account for them in perturbative QCD (pQCD) treatments. 
At a phenomenological level, the region of very small transverse momentum is most commonly described within a 
general transverse momentum dependent (TMD) parton model wherein the colliding hadrons are treated as collec- 
tions of nearly free point-like quark and gluon constituents. Within the parton model description, the TMD parton 
distribution functions (PDFs) and fragmentation functions are treated much like classical probability densities. They 
are conceptually very similar to the PDFs of the more familiar collinear parton model except that they describe the 
distribution of partons in terms of both longitudinal and transverse momentum components. 

In pQCD, a parton model description should be replaced by a factorization theorem. Demanding that factor- 
ization be consistent with real pQCD leads to rigid theoretical constraints which, in phenomenology, translate into 
first principles QCD predictions. When factorization theorems accommodate TMD PDFs and TMD fragmentation 
functions, they are usually called "TMD-factorization" or "/c^-factorization" theorems, and the constraints they im- 
pose are specific to cross sections that are sensitive to intrinsic parton transverse momentum. A derivation or proof 
of a TMD-factorization theorem for a particular process must show that the transversely differential cross section 
factorizes into a generalized product of a hard part (perturbatively calculable to fixed order in small coupling) , and 
a collection of well-defined non-perturbative TMD objects corresponding to TMD PDFs and TMD fragmentation 
functions for separate external hadrons. The non-perturbative TMD functions contain detailed information about the 
influence of intrinsic non-perturbative motion of bound state quarks and gluons, so measuring them and explaining 
them theoretically is currently the focus of many efforts to formulate a deeper understanding of hadronic bound states 
in terms of elementary QCD quark-gluon degrees of freedom. Because of their sensitivity to 3-dimensional intrinsic 
motion, TMD PDFs and fragmentation functions are also important in the study of the spin and orbital angular 
momentum composition of hadrons in terms of fundamental constituents. For lists of relevant references see, for 
example, Refs. [TJ[2]. In addition, for high energy collisions at the LHC, certain TMDs may be useful for studying the 
detailed properties of the Higgs boson H] . Generally, TMD-factorization theorems are necessary whenever highly 
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accurate calculations of transversely differential cross sections in the region small or zero transverse momentum are 
needed. This can include calculations used in the study of hadron structure as well as in searches for new physics. 

TMD-factorization theorems have been derived in pQCD with a high degree of rigor for a number of processes, 
including Drcll-Yan (DY) scattering, semi- inclusive deep inelastic scattering (SIDIS), and the production of back- 
to-back hadrons in e + e~ annihilation [5rtl2|. For those processes where TMD-factorization theorems are currently 
understood and known to exist, an important next step is to implement unified phcnomenological treatments, including 
global fits with QCD evolution and extractions of the non-perturbative TMD functions. The success or failure of these 
efforts will be a crucial test of the validity of small-coupling techniques in studies of hadronic structure and in new 
phenomenological regimes beyond what are treatable within ordinary collinear factorization. 

However, the focus of this paper is on processes where normal steps for deriving TMD-factorization in pQCD are 
now understood to fail. Recently, standard parton-model-based pictures of TMD-factorization have been found to 
conflict with pQCD for certain classes of high energy processes, most notably in high energy hadron-hadron collisions 
where a pair of hadrons or jets with large back-to-back transverse momentum is produced in the final state [13H19J . 
Interestingly, the conflict with TMD-factorization arises in kinematical regimes where common partonic intuition 
suggests that factorization should be very reliable. In particular, there is a hard scale, Q, set by the transverse 
momentum of the produced final state hadrons, which may be of the same order-of-magnitude as the center-of-mass 
energy. Therefore, the TMD-factorization breaking mechanisms discussed in Refs. [13H19] are distinct from other well- 
known kinematical complications with factorization such as those expected in the limit of very small Bjorken-x. They 
are instead due to an incompatibility between arguments for leading-power TMD-factorization and the non-Abelian 
gauge invariance of QCD, which persists even in the limit of a large hard scale. 

For processes and observables where factorization theorems are valid, Ward identities maintain factorization in the 
large Q limit by ensuring that any leading-power factorization breaking contributions that may appear term-by-term 
in perturbation theory cancel in the inclusive sum. After all cancellations have occurred, any remaining terms that 
violate factorization must be shown to be suppressed by powers of Aqcd/Q in order for factorization to be said to be 
valid at leading power. The details of these steps of a factorization proof constrain the specific form of factorization 
for the classic hard QCD processes. By contrast, the normally anticipated Ward identity cancellations fail in the 
scenarios discussed in Refs. [T3HTH], leaving leftover leading-power TMD-factorization breaking contributions. The 
non-cancellations were first noted in Refs. [131 E] for their ability to introduce interesting process dependence in 
TMD-functions, and they were later identified in Refs. [17l [18] as constituting a breakdown in the normal steps of 
a factorization proof. Still, it remained common to hypothesize that a more general form of TMD-factorization, 
called "generalized TMD-factorization" in Refs. jTHHIH], holds for these processes so long as the Wilson lines needed 
for gauge invariance in TMD definitions are allowed to have non-universaland potentially complicated structuresrj 
However, it was later found in [19] that TMD-factorization in the hadro-production of back-to-back hadrons fails even 
in this generalized sense. In other words, the TMD-factorization derivations cannot generally be made consistent with 
gauge invariance even when allowing for non-universal Wilson line topologies in the TMD definitions. Methods for 
addressing the resulting non-universality within tree-level or parton-model-like approaches have since been proposed 
in Refs. [20H23] . and methods tailored to the small Bjorken-x limit have been discussed in Ref. [24]l25j. 

A basic goal in the study of hadron structure is to establish a unified theoretical framework, rooted in elementary 
QCD quark and gluon concepts, for characterizing the intrinsic partonic correlations associated with QCD bound states 
and relating them to experimental observables. In processes where TMD-factorization theorems are known to be valid, 
the information about intrinsic hadron structure is contained in well-defined non-perturbative correlation functions 
like the TMD PDFs and the TMD fragmentation functions. The pattern that emerges is suggestive of a much more 
general picture of partonic interactions based on descriptions of quark-gluon properties for individual and separate 
external hadrons. With this broad picture as a foundation, the customary TMD classification schemes have usually 
assumed a form of generalized TMD-factorization for all types of interesting or relevant hard hadronic processes. The 
varieties of possible spin and angular behavior are then enumerated by first classifying the separate TMD functions 
according to their individual intrinsic properties, and then using them in a large set of both conjectured and derived 
TMD-factorization formulas. (We will elaborate in more detail on the meaning of "customary TMD classification 



schemes" in Sect. IV A ) 



This very general TMD-factorization picture has intuitive appeal because of the straightforward organizational 
method it implies and because it has a very direct connection to parton model intuition. In classifications of the 
possible spin and angular asymmetries of hard inclusive cross sections, it is almost always taken as an assumption. 
The only deviations that are usually allowed are those which incorporate non-universal normalization factors such as 
the overall sign reversal expected for the Sivers function in comparisons between the Drell-Yan process and SIDIS [55] 
or an overall color factor normalization for more complicated processes |16j . 



We will use the terms "gauge link" and "Wilson line" interchangeably throughout this paper. 



However, even the loosest forms of a TMD-factorization hypothesis impose strong constraints on the possible 
generalbehavior of transversely differential and spin dependent cross sections. Those constraints constitute valuable 
first-principles QCD-based predictions where TMD-factorization is expected to hold, but it has become common to 
also apply the standard classification scheme outside of the class of processes known to respect TMD-factorization, 
including in processes like those discussed in Refs. [T5HTU] . The purpose of this paper is to argue by way of example 
that, for these latter cases, the constraints imposed by a general TMD-factorization framework are too restrictivejj 
It will be shown that factorization breaking partonic correlations can produce unexpected patterns of spin and 
angular dependence in TMD cross sections that would otherwise be forbidden at leading power if a TMD-factorization 
hypothesis is adopted. Hence, the breakdown in TMD-factorization can modify the general landscape of leading-power 
spin and angular behavior in hadron-hadron scattering rather than simply giving process dependence to already known 
TMD functions. 

The relevance of TMD-factorization breaking extends beyond hadron or nuclear structure studies and is potentially 
important whenever perturbative QCD calculations are sensitive to the details of final states and are intended to 
have high precision point-by-point over a wide range of kinematics. For example, because of the detailed account 
of final state kinematics in TMD-factorization, it is a potentially valuable tool in the construction of Monte Carlo 
event generators [371 HSj- Also, it has recently even been found that factorization breaking arises in certain collinear 
cases [291 EQ] ■ In principle, it should be possible to connect the breakdown of TMD-factorization to the treatment of 
large higher order logarithms in collinear factorization and thus relate the two phenomena. 

Very generally, factorization theorems for inclusive processes rely on cancellations in the inclusive sums over final 
states. As such, the validity of any factorization theorem is placed in danger whenever extra final state constraints or 
conditions are imposed. A recent discussion in the context of collinear factorization for top-antitop pair production 
can be found in Ref. [3T] . For the TMD-factorization breaking scenarios of Refs. [T3HT§] , it is the specification of a 
small total transverse momentum for the final state back-to-back hadron or jet pair that is responsible for breaking 
TMD-factorization. Violations of standard factorization have been observed to have quite large phenomenological con- 
sequences, such as in measurements of hard diffraction in hadron-hadron collisions and in dijet projection, particularly 
when compared with measurements of hard deep inelastic diffraction [32H38] . 

The concept of a TMD PDF, and especially of an unintegrated gluon PDF, also appears in extensions of small 
coupling perturbation theory to the limit of small-x where other novel phenomena such as parton saturation become 
relevant. In this context, however, there are varying methods for using TMD gluon PDFs in calculations, and work 
is still needed to fully reconcile the different approaches with one another [24, 25, 39, 40 . 

A collection of noteworthy QCD-related phenomenological puzzles has developed over the past several decades. (See 
also Ref. [41.) This includes famously large and non-energy-suppressed transverse single spin asymmetries (SSAs) 
observed in hadron-hadron collision experiments at Argonne National Laboratory (ANL) [42l [43], Fermilab [111 
155] , CERN [43], Serpukhov [47], and Brookhaven National Laboratory (BNL) [4"5rl50] . For a recent review of the 
current experimental status of TMDs and spin asymmetries, see Ref. [5T]. Very early, large final state A-hyperon 
transverse polarizations were also observed in unpolarized hadron-nucleus collisions at Fermilab [52j . Similar transverse 
polarizations have not been detected in the inclusive e + e~ production of A-hyperons in the ALEPH experiment at 
LEP [53J, or in semi-inclusive deep inelastic scattering measurements by ZEUS [54] . More recently, an apparent sign 
discrepancy [55J has been identified in comparisons of Sivers asymmetries in SIDIS measurements [5BJ [57] with Sivers 
asymmetries extracted from h^h — > tt+X collisions at the BNL Relativistic Heavy Ion Collider (RHIC) [50] . Proposed 
explanations include a node in the Sivers function [55, [55], though recent fits disfavor this [59 . It is also possible that 
the apparent h?h — > 7T + X/SIDIS discrepancy is due to contributions beyond the Sivers effect [501 [5T]. At present, the 
full explanation remains unclear and more analysis is needed. At very high energies, interesting phenomena have also 
emerged, including the CMS ridge structure seen in proton-proton collisions at the LHC [62] ■ Explanations have been 
proposed in terms of Wilson line interactions [63] and the color glass condensate formalism J64H66] . In addition, a 
forward-backward asymmetry in the production of top-antitop quark pairs has been observed at the Tevatron J57J [55] , 
generating speculation about effects from physics beyond the standard model. 

Given the collection of observations in last few paragraphs, it is a natural stage now to investigate the possibility 
that TMD-factorization breaking effects lead to unexpected general patterns of behavior at leading power in the 
hard scale. This paper will demonstrate, by applying power counting to a specific example, that such an extra TMD- 
factorization breaking spin dependence can be induced by TMD-factorization breaking initial and final state soft gluon 
interactions. As our illustrative example, we will use the production of a final state large invariant-mass back-to-back 
photon-hadron pair in collisions of unpolarized protons. We will show that a single extra initial state soft gluon, of 
the type more commonly associated with gauge link effects, is sufficient to produce a correlation between the helicity 



2 The meaning of "generalized TMD-factorization," as it is used in this article, is much more general than what was called "generalized 
TMD-factorization" in Ref. [191 . The distinction will be explained in detail in Sect.lIVIof the main text. 
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FIG. 1. Diagram of the process in Eq. <TTj> in the coordinate system of frame-2. The final state hadronizing quark has transverse 
polarization Si and the final state prompt photon has helicity A3. 



of the final state photon and the transverse spin of the exiting final state hadron. We choose this specific example 
because it provides a relatively simple demonstration of an anomalous spin effect made possible by the breakdown 
of TMD-factorization while requiring only one extra soft/collinear gluon. However, the discussion is intended to 
illustrate the more general observation that spin correlations that would not ordinarily be expected to be large from 
the perspective a generalized TMD parton model can become leading-power in TMD-factorization breaking scenarios. 

The earlier examples of TMD-factorization breaking in Refs. [T51 - fPT)] focused on the appearance of anomalous color 
factors in low order graphs. In this paper, the TMD-factorization breaking will also involve an anomalous contraction 
of intrinsic transverse momentum components with Dirac matrices, and so it will appear to come from a different 
source. At the root of all these examples, however, is the breakdown in compatibility between non-Abelian gauge 
invariance cancellations and the factorizability of the transverse momentum dependent cross section at leading power 
in the region of small transverse momentum. 

The paper is organized as follows: In Sect. [TTJ we describe the process to be used throughout the paper to illustrate 
TMD-factorization breaking and the appearance of an extra TMD-factorization breaking spin asymmetry. Section In] 
also provides an overview of notation and conventions. In Sect. |III[ we give a brief summary of topics needed for 
discussions of TMD-factorization breaking, particularly in cases involving a TMD-gluon distribution in the initial state. 
In Sect. |IV| we precisely state the criteria for what we will call the maximally general form of a TMD-factorization 
hypothesis. Also in Sect. |IV[ we summarize the constraints that are placed on the spin and azimuthal dependence 
of differential cross sections by this maximally generalized TMD-factorization hypothesis. In Sect. [V] we present the 
detailed steps necessary to verify TMD-factorization at tree-level, and in Sect. |VI| we deal with spectator-spectator 
interactions. In Sect. | VII| we apply direct power counting to demonstrate a failure of TMD-factorization for the case 
of a single extra soft gluon attaching to the initial and final state, and we isolate the leftover TMD-factorization 
breaking contributions. We show in Sect. |VIH| that power counting implies an extra TMD-factorization violating spin 
correlation in the final state. In Sect. |LX]we summarize our main observations and make concluding remarks. 



II. SET UP AND NOTATION 



The Process 



For concreteness, we will focus on TMD-factorization breaking in the specific example of inclusive production of a 
back-to-back hadron-photon or jet-photon pair (see Fig. [I]) : 



Hadron 1 + Hadron 2 -» Jet(S r 4 ) + 7(A 3 ) + X . 



(1) 



+ Cross Terms 




FIG. 2. A parton-model level picture showing the 2 — >• 2 tree-level subprocesses qg — >• 97. 



The large transverse momentum of the final 7 and jet/hadron fixes the hard scale. The jet/hadron and the real 
photon are nearly back-to-back, so the cross section, which is differential in the total transverse momentum of the 
final state hadron-photon pair, is sensitive to the intrinsic transverse momentum of the constituent partons in the 
colliding hadrons. Hence, TMD-factorization is a natural candidate description. 

The outgoing quark in Fig. [Tlis meant to be representative of a final state hadron or jet with a transverse polarization 
04. 

Let Pi and P 2 be the momenta of the incoming colliding hadrons and let k 4 and k 3 be the momenta of the final 
state jet and real photon respectively. The total momentum of the final state pair is defined to be 



q = k 3 + k 4 , 



(2) 



and Q 2 = q 2 is the hard scale. The kinematical region of interest is characterized by final state transverse momenta 
of order the total center-of-mass energy of the colliding hadrons: 



(P + P 2 ) 2 - (P + k 3 ) 2 ~ (Pr + fc 4 ) 2 



0(Q 2 



(3) 



That is, the final state center-of-mass transverse momenta of k 3 and k 4 are large. 

The hard subprocess includes both qq — > 177 and gq — > 7*7 channels, but here we are most interested in the gluon- 
initiated subprocess since it is this that (as we will later show) most immediately yields a non-trivial extra spin 
asymmetry from TMD-factorization breaking effects. A parton model level depiction of the gq — > 707 hard scattering 
subprocess is shown in Fig. [2] the hard collision is the result of a quark inside Pi scattering off a gluon inside P2 . 
For a complete pQCD treatment, the qq — > 177 channel must also be accounted for, but an analysis of the gq — > 7(7 
subprocess is all that will be needed to obtain our main result. 

In the treatment of the hard subgraphs, we will often use the kinematical notation: 



§ =(h + k 2 y = (k 3 + k 4 y = q 2 = q a 

fc 4 - k 2 ) 2 = (fci 



t={k i -k 2 ) 2 = {k 1 -k 3 ) 2 =0(Q 2 ) 



(4) 
(5) 



We will also use a momentum variable defined as 



h = k 4 — k 2 = k\ — k 3 



(6) 



The proton mass is M p and the quark mass is m q . We are ultimately interested in spin correlations arising between 
the final state jet and prompt photon from factorization breaking effects, so we have labeled the jet/hadron and the 
photon with a spin vector, S4, and a helicity, A3, in Eq. (TlJ) and Fig. [I] 

A complete kinematical description of the cross section is closely analogous to the usual description of Drell-Yan 
scattering and e + e~ -annihilation into back-to-back jets. Our notation will therefore closely follow the standard Drell- 
Yan conventions. (See Ref. [12], chapter 13.2.) The modifications to the usual notation needed for the special case of 
Eq. (Til will be discussed below. 

First, we define the relevant coordinate reference frames: 




1. Pi, P* Center-of-Mass 

For categorizing the sizes of momentum components of the incoming partons, it will be most convenient to work 
in the center-of-mass of the colliding P1-P2 system, where Pi is moving in the extreme light-cone "+" direction and 
P2 is in the extreme "— " direction, and the total transverse momentum of the P1-P2 system is zero. We will refer to 
this as "frame-1." Momentum components in frame-1 will be labeled by an /1 subscript. The positive z-axis in this 
frame is along the direction of motion of Pi so that the four-momentum components of the incoming hadrons are 

p i./i = \Kt, > kS- . °* I « fe - 0. °*) . ( 7 ) 

(8) 
The total four-momentum of the final state jet-photon pair in frame-1 is 

1h = k 3J\ + fc 4, A = (<?/! > 1~h ' q *'/i ) • ( 9 ) 

In frame-1, therefore, small deviations of qtjj from zero are due to intrinsic transverse momentum in the colliding 
hadrons. The kinematical region of interest for studies of TMD-factorization and TMD-factorization breaking is where 
f\t,h i s °f or( ier A. (For the rest of this paper, A refers to a general hadronic mass scale like Aqcd-) 
The momentum fractions x\ and X2 arc defined as: 

k+ fi = Xl P+; k 2Jl = X2P2 • (10) 

The incoming partons k\ and k,2 are approximately collinear to their parent hadrons, so in frame-1 their components 
are of size 

k hfl ~ (Q, A 2 /Q, A*) , k 2Jl ~ (A 2 /Q, Q, At) . (11) 

In frame-1, the components of £3 and k^ are all large: 

fc3,A~(Q,Q,Qt), fc 4 , /l ^(Q,g,Q*)- (12) 

It will also be convenient to define unit four- vectors characterizing the extreme forward and backward directions in 
frame-1: 

ni, A s (1,0,00*, (13) 

n 2 , fl = (0,l,O t ) fl . (14) 

Most of the analyses of this paper will be performed in frame-1. 

At various stages, we will need only the frame-1 transverse components of a parton's momentum, expressed in 
four-vector form. Therefore, we define the four-vector notation: 

v t = (0,0,v t ). (15) 

The transverse components of the metric tensor in frame-1 may be conveniently expressed as 

9^ 2 = 9^ 2 - <A< 2 A - <,h<)h ■ ( 16 ) 

2. ks, ki Center-of-Mass 

We define "frame-2" to be the center-of-mass of the final state jet-photon system, and the corresponding momentum 
space coordinates will be labeled by subscripts /a . In this frame the orientation of the interaction planes is clearest 
and easiest to visualize. 

If the masses of the incoming hadrons (and the mass of the final state jet) are neglected, then in the center-of-mass 
of the jet-photon system 

<Z/ 2 = (Q,0), (17) 

fc 3 ,A = (<2/2,k /2 ), (18) 

fc 4 ,A = (Q/2,-k /2 ), (19) 

A,A = |Pi,/ 2 |(l,u li/2 ), (20) 

^2 )/2 = |P2,/ 2 |(l,u 2 , /2 ), (21) 



where here we have reverted to standard four-vector notation v = (v°, v 1 , v 2 , v 3 ). The energy of the jet-photon pair is 
Q, and u±j 2 , u 2 ,/ 2 are unit spatial three- vectors in the directions of Pi and P 2 respectively. In frame-2, the final state 
jet and real photon are exactly back-to-back, but the incoming hadrons P\ and P 2 are slightly away from back-to-back 
(see Fig. [I]) . The region of very small intrinsic transverse momentum corresponds to almost exactly back-to-back P\ 
and P 2 . The z-axis and x-axis can be covariantly fixed by defining the following four- vectors: 

Z n s (0»"i,/2- u 2,/ 2 ) X n s (0,uij 2 +u 2j2 ) 

|Ul,/ 2 -U 2 ,/ 2 | ' l U l,/ 2 + U 2,/ 2 l 

The z-axis bisects the angle between Pi,/ 2 and — P 2 ,/ 2 , and the z-axis is orthogonal to the z-axis, so that together 
they form a "Pi, P 2 -plane." Hence, frame-2 is analogous to the Collins-Soper frame [53] for Drell-Yan scattering, but 
with the lepton pair replaced by the jet-photon pair. The other plane, analogous to the lepton plane for Drell-Yan 
scattering, is formed by the jet-photon pair. We will call it the "£3, /cj-plane." An illustration of the process in the 
coordinate system of frame-2 is shown in Fig. [T] The polar angle of the jet axis with respect to the z-axis is 9f 2) and 
the azimuthal angle with respect to the Pi, P 2 -plane is <f>f 2 . 

3. &4 Boosted in the z-direction 

Finally, it will be useful in some instances to work in the ki, k% center-of-mass frame with the z-axis lying along the 
direction of fej, so that k± has its plus component boosted to order ~ Q and fc 3 has its minus component boosted to 
order ~ Q. We will call this "frame-3" and denote the components of vectors in this frame with an / 3 . In light-cone 
coordinates, 

*4,/3«(< /3 ,0,0) , (23) 

na, A = (0,1,0*)/,, (24) 

n 4> / 3 = (l,0,0t)/ 3 , (25) 

so that n 4 • k 3 = k^ , ~ Q. Frame-3 is simply related to frame-2 by a rotation. Most importantly for us, the transverse 
spatial components of Pi and P 2 in frame-3 are 

I (P lt ) /3 I « I (P 2t ) /3 I « |P li/2 | sm9 f2 = 0{Q) , (26) 

for small intrinsic transverse momentum and wide angle scattering (see Fig. fTl) . 

B. Final State Spin Dependence 

In later sections, we will be interested in the dependence on the polarization of the final state k^ in Fig. Ill and 
the representation of spin will follow the notation and conventions in appendix A of Ref. [12] . In particular, the spin 
vector Si of the final state quark in frame-3, where k^ is exactly in the z-direction, is 

= ( A 4fc4 j/3 sign(fc^ /3 ), m q b x tj3 ,m q b y thl \ A kl ih siga{kl h )) . (27) 

The A and b t j 3 are the components of the Bloch vector. The largest component of S4 is of order Q. The spin vector 
is normalized with a maximum of — S*| < m 2 ., and the spin sum is 



where kff ~ Q. We also define a vector 
so that 713 • &4 = k± t ~ Q, and 



E 



u(h, SiMh, Si) = (jfc 4 + m q ) (l +75— V (28) 

V m q J 



The polarization of the final state photon will be expressed in the usual way; the polarization four-vectors in frame-3 
are 

<-i = fo,-4=,— 5=><A , 4- 1 = (o,^=,-4=,o) . (29) 

A3 " X V V2 V2 ) h X3 — ' V V2 V2 ) h 

where A3 = ±1 refer to opposite helicities. 

For a review and a general treatment of spin see, e.g., Ref. [7D] and for typical conventions see [71] . 




FIG. 3. The structure of a general leading region. 



C. Leading Regions and the Classification of Subgraphs 

In perturbative QCD derivations of factorization in an arbitrary gauge, one must deal with extra soft and collinear 
gluons that are allowed in a general leading region analysis, but which do not factorize topologically graph-by-graph. 
Where factorization is known to be valid, the derivations show that contributions from these extra gluons either 
cancel in the sum over graphs, or factor into contributions corresponding to separate Wilson line operators, with each 
Wilson line operator belonging to a separately well-defined parton correlation function for each hadron. This standard 
factorization applies to sums over all graphs and relies on applications of the QCD Ward identity (gauge invariance) to 
disentangle the extra gluons that connect different subgraphs with one another. A factorization proof must show that, 
after the application of the gauge invariance arguments and Ward identity cancellations, any remaining leading-power 
factorization breaking terms cancel in the sum over final states for sufficiently inclusive observables. 

Standard leading region arguments [72] (see also chapter 13 of [73] and chapter 5 of [H]), show that graphs with the 
topology illustrated in Fig. [3] are the dominant ones at leading power for the process in Eq. (fTl). A generic Feynman 
graph may have arbitrarily many soft/collinear gluons (labeled lj, with j running from 1 to m) entering the hard 
subgraph Hl from the top gluon subgraph, $ s . A quark enters Hl from the lower quark subgraph, $ g , along with 
arbitrarily many extra longitudinally polarized gluons l[ (with i running from 1 to n). Analogous statements apply 
to the right side of the cut in Fig. [3j We define 



h = T,h 



(30) 



so that the total momentum entering the hard subgraphs, Hl and Hr, matches the routing of momentum in tree-level 
graphs like Fig. [2] There are also important contributions from gluons that directly connect the $> q and $> g subgraphs, 
but to keep the aiagram simple we have not shown them explicitly in Fig. [3] (see, however, Sect. VI I. The soft and 



collinear gluons may cross the final state cut, and this is included in the definitions of the $ ? and $ g bubbles. 

From here forward, the subscripts /i will be dropped, and it will be assumed to be implicit that we are working in 
frame-1 except when specified otherwise. In our notation, Lorentz indices denoted by Greek letters are intended to 
include all four components of a four-momentum, while indices denoted with i or j include only the frame-1 transverse 
spatial components. 

The differential cross section for an arbitrary number of extra soft/collinear gluons is expressed as 



Xi X 2 m,p,n,q 



d 4 k ± f d 4 k 



(2^)4 

d A U 



dH 

d% 



(2tt) 4 J (2tt) 



d% 



(27T) 4 

dHL 



d 4 l' 



(2tt) 4 J (2tt) 



d% 
{2-kY 



\M\ 2 x (mom. conserving 8— functions), (31) 




W(P 2 ; ft} )"'■"»' 



ff£(fci,fc3,fc4;{y;tt})»'' W1 ' 



£(Pi,*i;{y )«-■* 



FIG. 4. A generic contribution to the amplitude with the different pieces of Eq. (32 l labeled 



where M is the amplitude in Fig. HI and the integrals over the momenta entering or exiting the hard subgraphs are 
written explicitly. The inclusive sums and integral over final states in the upper and lower bubbles of Fig. [3] are 
represented by ^2 X Sx ■ The momentum conserving <5-functions fix k\ + fe = &3 + &4 — Q, with the definition of 
&2 in Eq. (30). For maximum generality, we leave the overall numerical normalization C unspecified in Eq. ( |31| . The 
essential requirement is only that the cross section is differential in q t — note that there is no integration over q t in 
Eq. (51). 

Inspection of the tree-level graphs (see, e.g., Sect. M) verifies that the power-law for M is M ~ Q°. The phase 



space integrals of Eq. (31 ) are Lorentz invariant and so do not contribute extra powers of Q to the cross section when 



boosting frames. In later sections we will analyze only the relative sizes, in powers of Q, of the factors that combine 
to form \M\ 2 . Therefore, the right-hand side of Eq. (31 1 is defined with a normalization factor Q 2 so that its overall 
power-law dependence is the same as for the amplitude, da ~ Q° (up to logarithmic factors). The overall Q 2 factor is 
to compensate for the two powers of 1/Q that come from the momentum conserving delta functions for k^ and k^. 
We are mainly interested in the dependence of the cross section on powers of Q, so to keep notation simple we will 
quote the power-behavior of any factor as a power of Q only, with any other factors of order a hadronic mass scale, 
necessary for maintaining correct units, left implicit. 

Since the cross section is order ~ Q°, then when a particular term in the perturbative expansion is of a higher 
power of Q than zero it is called "superleading." In a complete analysis including all graphs, gauge invariance ensures 
that all such terms cancel exactly against other superleading terms. 

In later sections we will find it simplest to work at the amplitude level and to consider the contribution to the cross 
section only in the last step. We will organize the analysis of amplitudes according to Fig. |4j A specific diagrammatic 
contribution can be broken into upper and lower subgraphs, U and C, and a hard part, Hi, as follows: 



M=U(P 2 ;{l j })^ 



Vifl ' ' ' fam X 

x u(k i: S A ) H L (k u *s, h; {l 3 } ; {J^iMi-AWi...*. x 

xg p[pi ---g plnPn C(P 1 ,k 1 ;{l i })^- p '«. (32) 



The Lorentz indices for gluons that connect U and C to the hard subgraph Hi are shown explicitly. In our notation, 
the final state quark wave function u(k±, $4) is written separately rather than being included in Hi because this will 
be convenient for power counting arguments in later sections. Hi and C also carry Dirac indices, though to maintain 
manageable notation, we do not show these explicitly. The subgraphs U, £ and Hi are defined to include the sum of 
all subgraphs corresponding to a given set of to external gluons for U, n external gluons for C, and to + n external 
gluons for Hi. However, this over counts the number of graphs by a factor of to! x n\. To compensate for this, 
therefore, we define the U subgraph to also include a factor of 1/to!, and the C subgraph to include a factor of 1/n!. 

All of our analysis will be done using Feynman gauge where the analytic properties of Feynman diagrams reflect 
relativistic causality and the main issues concerning TMD-factorization and TMD-factorization breaking are clearest. 
Lorentz indices will by denoted by /1, v, p, a, transverse components will be denoted by ij, and color indices will be 
denoted by a,/3,5,K. 

Actual derivations of factorization, where they are valid, involve a number of important points that will not be 
touched on in this paper, both because they are not directly relevant to the specific issues under consideration and 
because generalized TMD-factorization is anyway known to be violated for processes like Eq. (II]). Those other issues 
include, for example, rapidity divergences and other topological Wilson line issues that must be fully addressed in 
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situations where TMD-factorization is valid (see, for example, Refs. [74 ; -80 ), and their relationship to the evolution 
of separate, well-defined TMD PDFs |12j . Furthermore, in the integrations over intrinsic transverse momentum in 



Eq. (31), propagators in the hard part may go on-shell whenever any of the transverse momenta become too large. 
Therefore, in this paper it will be assumed that all such transverse momentum have upper cutoffs to prevent this, 
since the main concern is with transverse momenta of the order of hadronic mass scales or smaller. 

III. SUPERLEADING TERMS, GAUGE INVARIANCE, AND THE GLAUBER REGION 

A. Dominant Subdiagrams 

Recall that, in frame- 1, P± is highly boosted in the plus direction while P 2 is highly boosted in the minus direction. 
The graphs that individually contribute superleading terms to the cross section — i.e., with a power-law behavior 
higher than that of the cross section itself — are those in which all the components of IA are in the minus direction: 

U(P 2 ; {Ij})^-^ -> U(P 2 ; {lj})-'-'-'-'-. (33) 

The situation here is the same as was addressed in Ref. [81] in connection with ordinary collinear factorization and 
collincar gluon PDFs. Gauge invariance ensures that, in the sum over all diagrams, any superleading contributions 
to a physical cross section cancel. The way this is demonstrated relies on the same basic Ward identity arguments 
that are also needed to identify gauge link contributions and verify normal factorization theorems in processes where 
they apply [51152]. 

The observation from Ref. |81] that is most relevant to this paper is that, although the superleading powers of 



Q cancel after the sum over all graphs, contributions from subgraphs like Eq. (33) generally do not. (See, e.g., the 
third term in braces in Eq. (83) of Ref. |81|.) To understand this, recall that the minus component of polarization of 
the virtual fc2-gluon in hadron-2 is not exactly the same as its component along the direction of k 2 because k 2 also 
has order ~ A t transverse components. So a minus polarization is not exactly the same as a polarization along the 
direction of k 2l making the Ward identity cancellations more delicate. 



In spite of the non- vanishing leading-power contribution that arises from U subgraphs like Eq. (33), the cross section 
for the collinear case addressed in Ref. |81j does factorize, as expected, into the well-known and well-defined collincar 
gauge-invariant correlation functions. But even in that case the uncanceled leading-power terms from subdiagrams 



like Eq. ( 33 1 remain as contributions from the gauge link, and are therefore critical for maintaining consistency with 



gauge invariance at leading power. 

In later sections, we will encounter other complications from subgraphs like Eq. (33) in the treatment of generalized 
TMD-factorization breaking scenarios like Eq. (IT]) . 

B. Grammer-Yennie Method 

Factorization in covariant gauges is very conveniently addressed with the Grammar- Yennie approach [53] , which we 
briefly review here. In this method, each extra gluon's propagator numerator is split into two terms that we call the 
"G-term" and the "K-term." For a gluon with momentum I attaching U to the rest of the graph in Fig. [4J we write 
the decomposition as 



g ^V = K (iy n + G{iy» . (34) 



The G and K terms are defined as 



T7 P IV 

K{iy " = -^- (35) 

m ■ i 

G{IY " = g» " ■ - -±— . (36) 

n\ ■ I 

Note the explicit dependence of G(l) and K(l) on gluon four-momentum. The decomposition breaks the diagram into 
terms whose properties have a clear interpretation in terms of Ward identities and gauge invariance. The motivation 
for the separation is that it systematizes the identification of dominant and subdominant contributions to a subgraph, 
while automatically (and exactly) implementing any Ward identity cancellations that come from contracting the gluon 
momenta V?.-. with the rest of the graph in Fig. El Notice that in Eqs. (35 36) it is the exact l^ that appears in the 
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numerators rather than just l~ . After expressing the graphs in terms of G and K terms, it becomes clear how to 
apply valid approximations while at the same time maintaining gauge invariance. 

Equation ( 32 ) has been written with the metric tensors from the propagator numerators of extra gluons displayed 
explicitly. This is in preparation for later sections where we will implement the Grammer-Yennie decomposition on 
each extra gluon and examine th e powe r law behavior of all K and G combinations in each Feynman graph. 

If I ~ (A 2 /Q, Q, A t ) in Eq. (35 36) so that it is collincar, then power counting with the G and K terms is 
straightforward. For example, it is the g^ component of its propagator numerator in Eq. (32 1 that gives the 
dominant power, and the K^ term in Eq. (34) reproduces this exactly. As long as I ~ (A 2 /Q, Q, A t ), the other 

components of AT M M and G^ M are suppressed relative to the K + 



(34) reproduces this exactly. As long as I 

components. However, the simple power counting 



fails to apply directly if I is integrated into different momentum regions such as the Glauber region. 



C. The Glauber Region 



Many issues surrounding factorization and TMD factorization breaking are closely connected to the treatment of 
the Glauber region in the integrations over extra gluons in Fig. |4j so we briefly review the basic issue here. 

Consider, for example, a virtual gluon l^ attaching U and Hl in Fig. HI The transverse momentum may be of 
any size, though for addressing TMD-factorization we are especially interested in the contribution from the region 
where l t ~ A t . The integral over the l + component is then trapped by spectator poles at sizes of order ~ A 2 /Q. To 
recover factorization, I must either be strictly soft (all components of order A) or strictly collinear. However, in the 
integration over l~ there are generally leading regions where |/ + £~| << |^ 2 |- This is the Glauber region, where the 
Ward identities arguments that apply to ordinary soft or collinear gluons fail. To justify the steps of a factorization 
derivation, it must either be shown that it is possible to deform all contour integrals away from any Glauber regions 
so that they can be treated as properly soft or properly collinear, or the cross section must be sufficiently inclusive 
that any Glauber pole contributions cancel in the inclusive sum over final states. The potential for Glauber gluons 
to spoil factorization derivations has been noted long ago (e.g., Ref. [84]). 



IV. CONSTRAINTS FROM TMD-FACTORIZATION 

To ensure the clarity of later results, we must begin with a precise and detailed statement of what is meant when 
we speak of the usual TMD-factorization-based assumption as applied to Eq. ([l]) . The purpose of this section is to 
make a clear statement of what, for this paper, we will call the "maximally general" criteria necessary for a version 
of TMD-factorization to be said to hold. We will also summarize the constraints imposed by these assumptions and 
briefly discuss what it would mean for them to be violated. 

A. TMD-Factorization With Maximally General Criteria 

Our initial statement of the TMD-factorization hypothesis must be very general since the goal of later sections 
is to demonstrate that TMD-factorization fails for Eq. (fTl) even in a very loose sense. As reviewed in Sect. II C 



full derivation of TMD-factorization must account for the extra soft/collinear gluon connections of the type shown 
in Fig. [3j Thus, our statement of the maximally generalized TMD-factorization hypothesis is that, for N such extra 
gluons, in the sum over all graphs at a given order, the differential cross section in Eq. (31) must be expressible in 
the form 



da 



E / (2^5 ®9/r* (P*'M*»,<it - ku)) U Tr 



v 1 



u{k4,,Si) H l £ X3 (ki{xi 1 k 1 t),k 2 {x2 1 cit - k it ), /c 3 , fc 4 ) \ 

~i/Q 



ri, A 3 



x$ g / Pl [Pi;ki(xi,ku)) Hfi 3 [ki(xi,ki),k 2 {x2,q t - ki t ),/c3,fc 4 J u(ki,S4.\ 

,1/2 



i/Q 



0{{l/QY), (37) 



with a > so that error terms are power suppressed. 

We will now discuss the meaning of each factor in Eq. ( [37] ) . The basic structure is reminiscent of a TMD-factorization 
parton model, and it is a basic starting point for many phenomenological approaches. At this stage, the effective TMD 
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PDFs, &g/p 2 and ^ q /p 1 , need not necessarily correspond to gauge invariant operator definitions, within our maximally 
general criteria, and they may be highly process dependent. They may also carry any number of non-spacetime indices 
such as color (not shown explicitly). The polarization projections on Lorentz and Dirac indices for partons entering 
the hard subgraphs, H l £ 3 and H R ' 3 , are left as yet unspecified. Dirac indices will not be written explicitly. The M 
in front of the integration is intended to represent a possible collection of matrices in non-spacetime indices such as 



color, and it includes any overall numerical factors. The power law of Eq. (37) matches that of Eq. (31 ) for tree-level 
graphs like Fig. [2] With the horizontal braces beneath the equation, we have indicated the powers of Q for each 
factor, and we will use this notation throughout the paper to aid the reader in following the counting of powers of 
Qr] The way that the factors in Eq. (37) arise in a tree-level, parton model treatment will be reviewed in Sect.|v| 



The ordinary, bare-minimum kinematical parton model approximations needed to write down a factorized equation 
have been applied to initial state parton momenta in Eq. (371). That is, the minus component of ki is ignored outside 



its parent subgraph, ^ q /p 1 , and is integrated over within the definition of ^ q /p 1 (whatever that exact definition might 
be). Likewise, the plus component of ki is ignored outside & g /p 2 and is integrated over in the definition of & g /p 2 ■ The 
use of approximate parton momenta outside & q /p 1 and & g /p 2 is symbolized by the hats on ki, fc 2 and H L / R . Thus, 
the factorization applies after a mapping from exact to approximate k\ , 



fci^-fc 1 (xi I ki t ) = ^ 1 -(^,ki t ),fer(^,k lt ) J k lt (fe 1 + ,k lt )J , (38) 

is implemented outside $ 9 /p i; and a mapping from exact to approximate fc 2 , 

fc 2 -* h(x2,qt -k ti ) = (fej (fe^,k 2t ) ,%" (fe^,k 2t ) ,k 2t (fe^,k 2 t)) , (39) 



is implemented outside & g /p 2 ■ The components of the approximate hatted momentum variables, k^ and fc 2 , are written 
as functions only of fc+,k lt and k^ , k 2t respectively, as demanded by the minimal kinematical approximations. The 
only requirement needed in order for maximally generalized TMD factorization to be said to be valid is that the 
mapping from exact to approximate momenta be independent of k^ and fc 2 ; and that the approximate momenta 
obey the basic partonic four momentum conservation law, k\ + fc 2 = k\ + fc 2 = q. The precise definitions of the k\ and 
ki mappings may also depend on physical external hadron momenta, though this is not shown explicitly. The hard 



subgraphs in Eq. (37) have been replaced by their approximate versions, which use the hatted momentum variables: 

H L / R (ki,k2,k 3 ,k4,) -> B.ljr (fci(xi,ki t ), fc 2 (x 2 ,q t - k lt ),k 3 ,k 4 ) . (40) 

In a complete derivation of factorization, a basic step in the derivation would be to precisely formulate the exact- 
to-approximate mapping in Eqs. (38 39 ) such that a useful version of factorization is maintained to arbitrary order in 
small a s in the hard part. For our purposes, simply dropping the fcj" and and fc 2 components outside their respective 
parent subgraphs would be a sufficient prescription, but we have left the exact transformation laws unspecified in 



Eqs. (38 39) to ensure that the later demonstration of TMD- factorization breaking is completely general. Our only 



requirement for the mapping is that k x and k^ are neglected outside ® ip and $ fp ■ 

To summarize, we say that the maximally general TMD-factorization hypothesis is respected if, for N extra 
soft/collinear gluons in Fig. [31 

1.) ^ q /p 1 \Pi',ki{xi,'kit)) depends only on the momenta P\ and ki{x\,\L\t) and has the power-law 
behavior in frame-1 of an ordinary quark TMD PDF. It may be any matrix in non-spacetime indices, 
such as color, that contract with other factors in the factorization formula or with M. The Dirac 
polarization components of & q /p 1 may depend only on the target hadron momentum Pi and a target 
quark momentum k\. 

2.) & g /p 2 [P2] ki(xi,hit) j depends only on the momenta P 2 and fc 2 (x 2 ,k 2t = q t — ki) and has the 
power-law behavior in frame-1 of an ordinary gluon TMD PDF. It may be any matrix in non- 
spacetime indices that contract with the other factors in the factorization formula or with M. The 
frame-1 transverse Lorentz components, ij, may depend only on the total transverse momentum, k 2 t, 
of the incident gluon from hadron-2. 



3 



The way powers of large momentum ~ Q are partitioned between the hard factors and the effective TMD PDFs differs from some 
conventions — ours is chosen to simplify the power counting of later sections. 
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3.) H L / R (k 1 (xi,ki t ) , k 2 {x 2 ,'k2t) , k 3 , kA depends only on the momenta ki(xi, ku), &2(#2>k2t)j fe, and 
/C4, and has the power-law behavior of the tree- level hard subgraph. It may also be an arbitrary matrix 
in non-spacetime indices that contract with other factors in the factorization formula or with M. 



4.) The only exception to the first two criteria is that ^ q /p 1 ( -Pi; fci(a;i, kit) ) may include dependence on 
an auxiliary unit vector approximately equal to n 2l to account for possible gauge link operators, and 

<&g/p 2 ( P2; k 2 (x2, k 2 t) ) may include dependence on an auxiliary unit vector approximately equal to 



m. 

To be consistent with the maximally generally TMD-factorization assumption, the effective TMD PDFs and the hard 



parts in Eq. (37) may have complicated process dependence and may break gauge invariance as long as they satisfy 



these general conditions. 

Finally, note that the gluon TMD PDF in Eq. (37) only includes a sum over the transverse indices, i and j, for the 
gluon polarization, and not a sum over longitudinal indices. Within generalized TMD-factorization approaches, it is 
usually assumed from the outset that the basic polarization structure of the gluon subgraph is like that of an ordinary 
gluon number density, with only transverse components for a single on-shell gluon moving in the large minus direction 
contributing. Therefore, we include in the definition of maximally generalized TMD-factorization the assumption that 
effective gluon TMD PDFs always have two transverse Lorentz indices and depend on only one small transverse gluon 
momentum. 

Notice that the criteria enumerated above are less constraining than in what was referred to as "generalized TMD- 
factorization" in Ref. [19]. There it was required at least that separate gauge invariant (albeit process dependent) 
TMD PDFs be identifiable for generalized TMD-factorization to be said to hold. In this article, we require only a 



general separation into different transverse momentum dependent blocks in Eq. ( 37 ) with the properties listed above 



B. Standard Classification Scheme 

In this subsection, we review the conventional steps for extracting specific azimuthal angular or spin dependence 
in a TMD-factorization formalism based on the set of very general criteria laid out in the last subsection. The 
types of possible angular and spin dependence are attributed to special TMD PDFs that are defined with non- 
trivial polarization projections. The standard classification strategy, therefore, is to first enumerate the leading- 
power projections of Dirac structures in the effective unpolarized quark factor, ^ q /p 1 , and the leading projections on 



transverse Lorentz components in the effective unpolarized gluon TMD PDF $ fl /p 2 in Eq. (37). The result is a set of 
functions that are then identified with the various possible TMD PDFs. 

For the leading power unpolarized gluon TMD PDF, the only dependence on transverse momentum allowed by 



condition 2.) of Sect. IV A is from the total gluon transverse momentum k2t, and the only polarization dependence 
is from the transverse components i and j. Therefore, the leading power unpolarized gluon TMD PDF may be 
decomposed into the sum of a polarization-independent term and a polarization-dependent term: 

ij 

<& g /p 2 [P2;k 2 (x 2 ,q t -kit)' 



eff ^/P 2 (p 2 ;te,q*-kit))+ ^ + fp *S& (ft; *>(*»,% -k lt 



= P U *g/P 2 (ft! k 2 (x 2 ,q t - k lt )) + P B M(k 2t r $™ 2 (P2; k(x 2 , qt - k lt )) . (41) 

Here, §]!/p is the gluon TMD PDF for an unpolarized gluon and $Jm is a gluon Boer-Mulders [S5J [SBJ function. In 
the last line, we have defined projections onto transverse indices: 

^=~9?, (42) 



PBM(k*)« = ( M* + fLg* ) . (43) 



Ui jj J-3-L.2 

^2t^2t 1 tit ft : 
An auxiliary gauge-link vector n\ adds no other structure since it has no transverse components 
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The unpolarized quark TMD PDF, <& q /p 1 , is decomposed into a complete basis of Dirac matrices, 



*,/*, [PL'Mxuku)) =S + lb P + 1 ^+ 1 ^ b A^ + -a^T^ ~Q ; 



(44) 



and we keep only the leading power components. Since & q /p 1 ( P\ ', k\ {x\ , kit) ) depends only on Pi and fci(:ri,k lt ), 

which in frame- 1 have large (order Q) components in the plus direction, then the leading terms in Eq. (44) are those 
with a fj, = + contravariant index. The TMD PDFs are therefore obtained from Dirac traces that project those large 
components. For example, the largest component of the vector term, j^V^ , is 

7 -V + ~Q. (45) 

The ordinary azimuthally symmetric and unpolarized TMD PDF is therefore identified with the projection 



*[/2 (Pi-MauM)) =V+ = ^Tr [ 7 + *,/ft (Pi;*i(a;i,ki) 



Q. 



(46) 



In this expression, we have adopted the notation of Ref. [57] by including a [ 7 + ] to indicate the specific Dirac projection. 

U+l 
The only dependence on external momenta is through Pi and ki, so the largest Dirac components of $ , p are 

proportional to P ~ $ 1 ~ Q. The trace in Eq. (46) then takes the general form, 



$ [ 7 /p! (Pi;fci(^i,ki)) ex ^Tr [ 7 +A] ~ Q. 



(47) 



An auxiliary gauge-link vector i/l 2 introduces no further leading-power projections because it only leads to additional 
factors of 7 + in Eq. ((47 1 . 



For a projection with a generic Dirac structure T q , the corresponding quark TMD PDF is 



*J^ (Pi;*i(a;i,ki)) = ^Tv[r q <P q/Pl (^P 1 ;k 1 (x 1 ,k 1 ) 



Q. 



(48) 



The decomposition in Eq. (44) spans the set of possible independent Dirac structures for a quark TMD PDF, and 



each T q is associated with a characteristic type of angular or spin dependence in the cross section. In the TMD- 
factorization-based classification scheme, the types of leading-power quark TMD PDFs are categorized according to 
the possible Dirac projections T q that contribute at leading power. For example, the projection, 



iW 7 + 



(49) 



reproduces the ordinary azimuthally symmetric and unpolarized TMD PDF in Eq. (46). Another well-known example 
is the quark Boer-Mulders [85 function, obtained from the projection, 



r g -> 7 + y 7 5 . 



(50) 



If hadrons 1 and 2 are also polarized then there is a large collection of additional leading-power TMD PDF structures 
involving hadron spins. These have been thoroughly classified at least to leading twist [571155] , 

So the maximally general statement of the TMD-factorization conjecture of Sect. IV A when specialized to a 
particular set of projections, T g and T q , is 



daPs^ 



E 



^$[ r /Y i 2 (p 2; fc 2(X2 , qt _ klf 



(2ny 



$ 



[r„] 

9/ Pi 



Pi;fei(xi,kit) 



xPr s ( qi -k lt )«Tr 



l(k4,S 4 ) H l £ Xa (k 1 (x 1 ,k lt ),k 2 (x2,q.t~'kit),k 3 , 



k 4 ] x 



-,1/2 



VO 



x Pr, H R 



j,h 



[k 1 (xi,k lt ),k 2 (x 2 ,qt -k 14 ),fc 3 ,fc 4 j u(fc 4 , 5 4 ) 



(51) 



I/O 



This is obtained by using Eqs. (41144148) in Eq. (37). The Pr g on the last line of Eq. (51 1 is the relevant Dirac matrix 
from Eq. (44) corresponding to projection T q . (In the unpolarized azimuthally symmetric case, Pr c — > 7 — •) The 
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r g = {U,BM} labels which gluon polarization projection is taken from Eq. (41). The superscript [r ? ,r g ] on da is to 



indicate that Eq. (51) is the contribution to the cross section from that particular combination of T q and T g target 



•J r * 



iP. 



quark and gluon polarization projections. 

Throughout this section, we have used the language of a generalized parton model by referring to 3> L± /p and $ l , p 
as though they were ordinary TMD PDFs. However, in general pQCD diagrams like Fig. [3j the $'s are merely labels 
for the target subgraphs corresponding to hadrons 1 and 2, each of which may be linked to the rest of the process via 
arbitrarily many soft and collinear gluon exchanges. Strictly speaking, TMD PDFs acquire a well-defined meaning 
in pQCD only after a factorization derivation has shown how those extra soft and collinear gluons separate into 
independently defined factors in the sum over all graphs. 

Because of its very general form, there is a strong temptation to apply the factorized classification scheme of 
this subsection very broadly, even in processes like Eq. (fTl) where generalized TMD-factorization formally fails |19j . 
But even the maximally general version of the TMD-factorization conjecture from this section imposes significant 
constraints on the possible general behavior of the cross section at leading power. For example, the power counting 
logic used to identify the leading TMD functions in Eq. (41) and Eq. (44) relies crucially on the assumptions that, 
a.) [*„ 



'g/P-2 



depends at leading power only on transverse Lorentz inc 



ices and the momenta Po, k and k? f and, 



b.) ^? q /p 1 depends only on momenta Pi, k{~ and ki t . When generalized TMD-factorization is broken, we must in 
principle account for the possibility that dependence on extra external momenta leaks into the & g /p 2 and & q /p 1 from 
other subgraphs via the extra soft and collinear gluon exchanges in Fig. [3j However, in confronting this possibility we 
are no longer justified in assuming from the outset that, for example, the [i = + components in Eq. (H4j) are the only 
dominant ones. Any Dirac structure projected from Eq. (44) must be viewed as a potential candidate leading-power 
effect. Analogously, the decomposition of [& g /p 2 Y : ' may include terms beyond those of Eq. (41 ), involving polarizations 



induced by internal intrinsic transverse momenta other than k2t, and polarization configurations like Eq. (33) must 
also be accounted for. We will discuss the implications of this in more detail in the next subsection. 

For the sake of clarity, our discussion has so far focused on the TMD PDFs, though the same observations apply 
to a treatment of final state hadronization, i.e. to the treatment of a quark fragmentation or jet function. The 
use of TMD-factorization, in the very general form summarized in this section, has become standard for classifying 
leading-power TMD effects. It is a valid method, of course, in processes where a form of TMD-factorization holds or 
is expected to hold (such as in SIDIS or Drell-Yan). For TMD-factorization breaking processes, however, we will find 
that the structure in Eq. (|51|) turns out to be overly restrictive. 



C. Parity, TMD PDFs and TMD-Factorization Breaking 



A TMD-factorization conjecture imposes even stronger constraints on the qualitative structure of cross section when 
considered in combination with discrete symmetries of QCD. Discrete symmetry arguments have long been useful for 
constraining the properties of parton correlation functions in pQCD. It is now known, however, that such arguments 
are more subtle in processes that involve TMD-factorization than in similar collinear cases, due to the non-trivial 
role of gauge links in TMD-factorization. Early on, for example, TP invariance was applied in Ref. [89 to argue 
that the Sivers function should vanish in hard processes. That derivation, however, neglected effects from gauge links 
in the definitions of TMD PDFs. As the role of gauge links in the derivation of TMD-factorization became better 
understood, it was realized that TP invariance implies not that the Sivers function vanishes, but rather that it reverses 
sign between the Drell-Yan process and SIDIS [351 EUl HI] > thus providing a non-trivial prediction for a specific type 
of non-universality 

Such observations provide motivation to also examine how arguments based on discrete symmetries arc affected 
when TMD-factorization breaks down altogether. As an example, let us return to the quark TMD PDF's Dirac 
decomposition in Eq. ( 44 ) . In addition to Eq. ( 49 ) , naive power counting would imply another leading-power Dirac 



structure formed by the axial vector term 

757" A+ -Q. 
So one should also include in the categorization of leading projections, 

r 9 -> 757 + ■ 



Therefore, by analogy with Eq. (46), a TMD function defined as 

1. 



47ft (^i;*i(*i.ki)) 



-Tr 



757 ®q/Pi (-Pi;fci(#i,ki, 



Q 



(52) 



(53) 



(54) 
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should be included within the classification of leading-power quark TMD PDFs. The interpretation of this TMD PDF 
in a generalized parton model/TMD-factorization approach would be that it describes a distribution of quarks with 
an intrinsic helicity polarization inside an unpolarized hadron. Such TMD PDFs are forbidden by parity invariance 
for an unpolarized parent hadron for any TMD PDF definitions that satisfy the basic requirements of the generalized 



TMD-factorization hypothesis enumerated in Sect. IV A In Feynman graph calculations, this appears in the inability 
to construct a pseudo-scalar using only the four-momenta from hadron- 1. At least four different four-momentum 
vectors are needed to give a non- vanishing trace. Specifically, an evaluation of Eq. (54) always leads to a Dirac trace 
of the form, 



$ 



[757 J / 
q/Pi \ 



P t ;ki(xi,l<-i)) ex -Tr 



757 



-ktfi 



0. 



(55) 



There are not enough four momentum vectors involved to give a non-vanishing trace. A gauge link simply introduces 
more ~ rii vectors, which only project extra factors of 7 + in the Dirac trace: 



$ 



[757+] / 

q/Pi \ 



Pi;fei(a;i,ki)) ex -Tr 



ISl+htf^ 



^Tr 

4 



757 



'ktfi 



7 



0. 



(56) 



So even including gauge links does not allow for projections of this type within a TMD-factorization conjecture. TMD 
PDFs like Eq. (55), and any effects that would follow from them, are constrained to vanish at leading power within 
the minimal TMD-factorization criteria of Sect. |IV A"| and |TVB| For this reason, they might be referred to as "naively 
parity violating." 

To further investigate the consequences of discrete symmetries in a generalized TMD-factorization framework, we 
will focus in later sections on the final state parton polarizations — the helicity A3 of the prompt photon and the 
the transverse spin S^x of the outgoing quark in Fig. [T] If the initial state partons have no helicity then the hard 
subprocess 



■,(Sa 



+ 7 



(A3) 



(57) 



violates the helicity conservation of massless QCD. For the S4 j_ — A3 correlation to be possible, a helicity must be 
carried by the q or g inside its parent unpolarized hadron before the hard collision. In a generalized TMD factorization 
framework, this would imply a type of naively parity violating TMD PDF. 

After showing that the criteria of Sect. IV A fail to hold in a more detailed treatment of the extra soft/collinear 
gluons of Fig. [3) we will ultimately demonstrate in Sect. VIII that final state S^x-^3 correlations are actually leading 
power in spite of the apparent helicity non-conservation in the naive 2 — > 2 partonic description in Eq. ( |57[ ) . 

If one abandons the maximally general TMD-factorization criteria of Sect. |IV A| and allows the effective TMD 
PDFs to acquire dependence on additional extra external momenta, then any Dirac structure projected from Eq. (44) 



is a potential candidate for a leading-power quark TMD effect, rather than only those with a "+" Lorentz component. 
For example, the scalar term, 5, in Eq. (44) is at most of order a hadronic mass within the TMD-factorization 
framework of Sect. |IV A| and so would be counted as subleading. But if, because of the extra gluon connections in 
Fig " 



§ fq/Pi 



Pi; k\(xi, ki) J were allowed to depend on other external momenta, 

*q/Pi (Pi;fcl(»l,kl)J — > * q /Pi yPi;ki(xi,ki),q,P 2 ,k 2 . 



(58) 



then S could acquire dependence on large scales like Pi • Pi and q ■ kit- Then the usual power counting arguments for 
Qq/p-L become invalid. 



The pseudo-scalar term, 75P, in Eq 
invariance in the TMD-factorization framework 



(44) is another example that is both subleading and forbidden by parity 
It will be the main example that we will use in Sect. |VHI| to 
demonstrate a leading-power TMD-factorization breaking final state spin correlation. The projection is simply 



and the substitution in Eq. (51 ) is 



r 9 ^75 



r a ->76 

= 75 ■ 



(59) 



(60) 



If one works within the TMD-factorization conjecture by following the classification scheme of Sect. IV B then one 
would define a TMD PDF corresponding to Eq. (48) using the projection in Eq. (59): 



$ 



[75] 
q/Pi 



Pi;fci(xi,ki 



1. 



-Tr 



l^q/Pi (Pi;fci(a;i,ki)J cxm g Tr[7 5 ] =0 



(61) 
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So, the possibility of any spin correlation effects arising from such projections would be discounted in the standard 
TMD-factorization classification scheme. 



We will return again to the projection in Eq. (591 in Sect. VIII where we will confront TMD-factorization breaking 
effects. 



D. Maximally General TMD-factorization at the Amplitude Level 



The conditions enumerated in Sect. IV A are sufficiently weak that it becomes straightforward to separate, at the 
amplitude level, contributions that are definitely consistent with the maximally general TMD-factorization criteria 
from those which may violate TMD-factorization if left uncanceled. This will be useful for later sections where much 
of the analysis is performed at the amplitude level. 



Assume that for a fixed number m and n of extra gluons in Eq. ( 32 ) one is able to sum a set of Feynman diagrams 
to obtain a contribution to the amplitude of the form 



M 



■■M J2 U eff(P2-,k 2 ,{l 2 , 



,,i m }y u(k 4 ,s 4 ) Hl^mMM 1 • 3 £eff(Pi;h,{i[,...,i' n }) 



(62) 



i/Q 



Here we have changed variables from l\ to k 2 = ^7J ■ lj- The factor TV can represent any matrix of non-spacctime 
indices such as color. The hard subgraph H^ depends, as in Sect. |IV A[ only on the momenta ki, k 2 , fc 3 and k 4 . 
Assume further that the effective U e jf and C e ff factors take the forms, 



Z4//(P 2 ; k 2 , {l 2 , • • ■ , Ira})' = T U (Pl\ k 2 , {h, ■■■, lm}Y 



V u (k 2t Y 



and 



£eff(P 2 ; k 2 , {l[, . • . , G) = FciPx; hi, {l[, ..., 0) p; " p " Pc{PiMt. 



■ci, 



p'v-p'-n. 



(63) 



(64) 



By substituting Eq. (62) into Eq, (|31[) and applying the minimal kinematical approximations of Eqs. (38M0), one 



immediately recovers the structure in Eq. ( 37 ) 



Note that Tu and Tc do not share any ^-momentum arguments. Since there is no direct dependence on {Z 2 , . . . l m } 
outside U e ff , the part of Eq. ( 63 ) labeled by Tu , which does not depend on transverse polarization components (labeled 
"i"), may have dependence on any of these extra gluon momenta and still maintain consistency with the maximally 
general TMD-factorization criteria of Sect. |IV A") In the squared amplitude, these momenta are simply integrated over 
inside & g /p 2 - However, the polarization-dependent part, labeled by VuQ&2t) % ■> couples to the transverse "i" component 
of the hard part in Eq. ( 62 ) , so it is allowed to depend only on the total transverse momentum, k 2t , coming out of the 
gluon subgraph. Without such a requirement, the criteria from Sect. |IVA| can be violated, and the decomposition of 
the effective gluon TMD PDF into polarization dependent functions in terms of k 2t , as in Eq. (41), would in general 
be incomplete. The integrations over \ 2 t ■ ■ ■ could not be performed independently of the contractions in the hard 
part. 

Similarly, there is no direct dependence on {l[,... l' n } outside C e ff, so the non-polarization dependent part of 
Eq. (64), labeled by J 7 ^, may depend on any of these extra gluon momenta. In the squared amplitude, they will be 



integrated over inside the definition of ^ q /p 1 ■ However, the projection onto Dirac components, labeled by Vc(Pi, kit), 
can only depend on Pi and the total momentum k\ of the quark coming out of hadron-1. 

A potential violation of the criteria of Sect. |IV A| would arise, for example, from a contribution in which one is 
forced to write Eq. ( 63 1 in the form, 



U ef f(P 2 ; k 2 ,{l 2 , ..., l m }) 1 = T U (P 2 ; k 2 , {l 2 , . . . , l m }) 



Ml'" Mr, 



^w(l2t,k 2 t,...)^..x 



(65) 



i.e., where the polarization projection Vu depends separately on the transverse components of individual extra gluon 
transverse momentum, rather than only on their sum. Recall that it is only h 2t that enters the hard scattering. 
In Eq. ( 65 ) the transverse gluon momentum k 2 * that enters the hard subgraph may differ from the transverse gluon 
momentum \ 2t that induces a polarization dependence. If such terms fail to cancel, then they contribute to a breakdown 
of TMD-factorization, even under the maximally general criteria of Sect. |IV A| 



V. ANALYSIS OF ONE GLUON 



Before addressing the case of multiple gluon interactions, we will use this section to give a detailed demonstration of 
how Ward identity arguments apply to the simplest case of just one gluon attaching hadron-2 to the rest of the graph 
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(a) 



(b) 



FIG. 5. The single gluon parton-model-level contributions to Fig. HI 



(Fig. [5]). In the absence of extra gluons, there are few enough complications that the problems with TMD-factorization 
will not appear — extra soft and collinear gluons of the type normally associated with gauge links are needed for the 
violation of TMD-factorization to be apparent. Therefore, the results of this section will be found to be consistent 
with the conditions for maximally general TMD-factorization from Sect. |IV| Nevertheless, the steps will help establish 
the basic framework needed for dealing more generally with factorization issues in Sect. |VH| For simplicity, in this 
section we will restrict consideration to the case where the final states in Fig. IT] are totally unpolarized. 



A. Separation into Subgraphs 

We begin by applying the analysis of Sect. [In] to the special case of Fig. [5] The separation into subgraphs according 
to Eq. p2| is 



M {lg) =U(P 2 ;k 2 ) K ' ti J^SW H L (ki,k 2 ,k 3 ,k 4 r<^ £(P l5 fc!) 



(66) 



The Lorentz index for the gluon is /i, that for the photon is a, and n is the color-octet index for the gluon. Repeated 
color indices are summed over. We will drop the explicit a index on Mi lg \. 

In frame-1, the powers of Q for the components oili^ are 



U~ 



Q, U l ~ Q°, 



ir 



l/Q. 



(67) 



As usual, superscripts u i,j" denote transverse Lorentz components in frame-1. The final state momenta, k% and ki, 
are at wide angles so all components of H£ are of comparable size: 



H-, 



m 



m 



l/Q. 



(68) 



For this section, the Lorentz index on Hl and IA always refers to the gluon unless otherwise specified. 

The final state quark wave function behaves as u(k±) ~ Q 1 ' 2 . The dominant contribution from C comes from 
leading-power projections onto the Dirac components of $> q /p 1 , as in Eq. (46 1, which are at most of order Q. So the 
dominant component of C also has a power-law behavior C ~ Q 1 ' 2 . 

Since we are mainly interested the contraction of hi with the rest of the graph, let us introduce the notation, 



R K '^ = u(ki) H L (fc 1 ,fc 2 ,fc 3 ,fc4) K; ^ C(P r , fci) , 

-Q 1 / 2 ~1/Q -Q 1 / 2 

so that 

A/ (ls) =W(P 2 ;fc 2 r iM ' j^iJ s; '" • 

Then all the components of R K '^ (T have the power-law behavior, 

RK ;-a _ R n;+a _ RK ;io _ qO _ 



(69) 

(70) 
(71) 
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and the largest power-law in the full amplitude from the tree-level graphs of Fig. [5] is from the contraction 

M ™g) ~ U K ~R K ' +,T ~ Q . (72) 

This is a contribution from the X-gluon term in the Grammer-Yennie decomposition. It is a superleading contribution 
of the type discussed in Sect. |HI| 

We next implement the Grammer-Yennie separation in Eqs. (34f36) in Eq. (70 1 by writing 



M {lg) = M, 



(is) 



M, 



K 

(is) 



with 



and, 



Mf Xg) = U(P 2 ; k 2 r^'G(k 2 )^,R K ^ 
-U(P 2 ;k 2 ) K ^K(k 2 ),^R K ^ 



M, 



(is) 



G-term , 
K-term . 



K(k 2 ) 



n i k 2 

m ■ k 2 



G{k 2 y " = g» 



l i k 2 



(73) 

(74) 
(75) 

(76) 
(77) 



rii ■ k 2 

The if -term gives superleading as well as subleading contributions in the amplitude. The G-term, however, has been 
deliberately constructed so that its "H — " component is removed. The highest power-law contributed by the G-term 
when contracted with U is therefore 



7 jK',i T3K\lO 



U K <- R K:ia {k 2t /k^) ~ Q°. 
~ Q ~ l/Q 



(78) 



These are the leading-power G-term contributions to the amplitude. The remaining G-term contributions, in the 
contraction of U with R, are power suppressed relative to the leading terms: 



i/Q, 



(79) 



in accordance with the power-laws in Eqs. (67 1 and (71). The explicit expression for R K '^ a , from the two graphs in 
Fig. § is 



/« - ig s e q u{k A ) 1 {h h)2 +iQ + {h + h)2 _ + , Q J» C{Pi, ki) 



(80) 



B. Cancellation of Superleading Terms 



noting that Eq. (|75|) can be written as 



The superleading K-term in Eq. ( 70 ) can be seen to vanish in the sum over graphs in the full amplitude by first 

_u{p 2 -h 2 y- 



U(P 2 ; k 2 )^ K(k 2 ) ll > ll R Kl '" T = U{P 2] k 2 ) 



■ni ^k 2 ^ 
ni ■ k 2 



R K 



X k^R^" ~Q, 



(81) 



so that attention may be focused solely on k 2tl R K ' tlc ' . (The underbraces in Eq. (81 ) only denote the maximum powers 
of Q graph-by- graph. Since they are superleading, the terms of order Q have to cancel in the sum over all graphs.) 
In order to argue that all the if-gluon contributions in Eq. (81) are power suppressed, it must be shown that both 
the order ~ Q (superleading) and the order ~ Q° (leading) contributions to the k 2u R K ' ,fi(T subgraph cancel in the sum 
over graphs when contracted with hi(P 2 ; k 2 ) K '^ . Including both terms from Eq. (80), fc 2p i? K;A " J becomes 

k 2il R K ^ a = ig s e q u(h +k 2 - fc 3 ) 

h(h-h+m q )rt K , Terml 

Term 2 



(ki — k.3) 2 — m 2 , + iO 

r{h + h+m q )ht R \ 

(ki + k 2 ) 2 -m 2 + i0 
xC(Pi;ki). 



(82) 
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From here forward, the steps for dealing with Fig.[5]are similar to standard arguments for a Ward identity cancellation. 
For Term 1 in braces in Eq. (82), we apply the Feynman identity replacement, 



while in Term 2 we use 



When the 



r2 — (rl + r2 P3 m <j) \rl $3 rn q) ; 



¥•2 ~ \r"X "I" r-1 m q) Vri m <?) ■ 



(83) 



(84) 



m q ) in Eq. (83) acts to the left on u(ki + k 2 — fc 3 ), in Term 1 of Eq. (82), the resulting 



contribution vanishes exactly by the action of the Dirac equation on the outgoing quark wavefunction. 

When the — ($i — m q ) in Eq. (84) acts to the right on C(Pi;ki), in Term 2 of Eq. (83), the resulting contribution 
is suppressed by two powers of Q. (It would be exactly zero if the target quark were taken to be exactly on-shell.) 
To see this, let us rewrite Term 2 from Eq. (82), but with the ki -propagator inside £(Pi;fci) explicitly displayed: 



ig s e q u(ki + k 2 ~ fc 3 ) 



r(h + h 


+ m q )$ 2 t R 


(*i + h) 2 


— m q + iO 



C(P 1 ;k l ) = 



ig s e q u(ki + k 2 - k 3 ) 



ri.h+h+m q y 
(k 1 + k 2 ) 2 -r l 



iO 



($! +m q ) 



"2l,2 



k\ — to 2 + iO 



Q . (85) 



i/Q 



In the last line, the "• • • " symbolizes the other factors that make up £(-Pi; fci), and these are order ~ Q~ x / 2 because 
of the overall C{P\\ k\) ~ Q 1 ^ 2 power-law. Note that in the second line of Eq. (85), the ft 2 has been moved outside 
the hard propagator and is included in the parentheses with the fci-propagator. 



When the Feynman identity substitution is made for $ 2 , the — ($ l — m q ) from the second term of Eq. (84) gives 



ig s e q u(ki + k 2 - k 3 ) 



,1/2 



r{h + h+m q )t K 
(fci + fc 2 ) 2 -m^+iO 



x [-(^-m,)] 



k\ - m 2 + iO 



VQ 



(86) 



i/Q 



Q° 



So the factor in parentheses looses two powers of Q relative to Eq. (85), and the contribution to k 2ii R K '^ a from the 
second term in Eq. (84) has a power law ~ 1/Q. That it is suppressed by two powers of Q instead of only one is 
crucial because we will ultimately multiply it with the superleading U~ in Eq. (81 1 . We will need this result again in 



Sect. VII when we take into account effects from an additional soft/collinear gluon radiated from hadron-2. 

The remaining if -term contributions are from the second term in Eq. (|83|) and the first term in Eq. ( 84 ) in the 



Feynman identity substitutions. These exactly cancel against each other when we use quark propagator denominator 
cancellations like 



(fci + k 2 y - 



m q ) 



to 2 . + iO 



(h 



-m q ) = 1. 



(87) 



The exact cancellation is not put in danger by the «0 because (fci + k 2 ) 2 is constrained by kinematics to be always of 
order Q 2 . 

Thus, there is no leading-power if-gluon contribution in the one-gluon case, and the superleading contributions 
cancel as expected. The only remaining contribution is from the G-gluon in Eq. (74). From Eq. (78), it immediately 
follows that the leading-power G-gluon contributions to Mn g ) involve only the transverse /^-components of i? K;/i<T . 
The plus component is already entirely accounted for above in the treatment of the if-term and is, by construction, 
exactly removed in the G-term. Because of the suppression in Eq. ( 79 1 , the minus /i-component of R K, '<i ia gives a 



contribution suppressed by a power of Q relative to the leading power. Therefore, the Lorentz components in the 
G-term contraction may be restricted to the frame-1 transverse /x-components. That is, we may replace 



G(k 2 y 



G(k 2 y 



k. 
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1%\ ■ k 2 



to leading power, which is in agreement with the normal gluon PDF vertex with no gauge link contribution (see, e.g., 
Eq. (45) of Ref. [5T| and Fig. 3.4 of Ref. [B]). The lowest-order uncanceled contribution to the amplitude is therefore 
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consistent with typical expectations for the parton model gluon density: 



M (ls) = z«? s e 9 > U(P 2 ;k 2 );, 5 "* - 



rii • fc2 



w(fc 4 



J I (fci - fc 3 ) 2 - m2 + iO (fci + fe) 2 - to 2 + zO J l *' lj 



'.9. 



,e,y)Z/(P2;*2)*, G(k 2 r' i ufa) H h o(k u k 3 ,k 4 ,k 2 ) K ^ £(Pi;fc x ) ~ Q° , 



i/Q 



(89) 



where we have used Eqs. (88 80) in Eq. J74[ ). In the second line, we have defined H^oiki, k 3 , k 4 , k 2 ) K '^' 7 as the sum 
of the lowest order contributions to the hard parts in Fig. [5] 



Hho(ki,k 3 ,k4,k 2 ] 



K,;fi(T 



(hi - h) 2 -m 2 + i0 (fci +fc 2 ) 2 -m 2 + i0 



(90) 



Checking Eqs. ( 89|90 ) with Sect. IVD confirms that the amplitude is now in a form consistent with the TMD- 
factorization conjecture, Eqs. ( 62]|64 ), in the maximally general form. Tallying the powers of Q indicated with braces 
underneath Eq. (89 1 verifies that the combined powers give a leading contribution to the amplitude. 



C. Tree Level TMD-Factorization 

Now nothing obstructs the immediate recovery of the normal TMD-factorization formula characteristic of parton 
model expectations. Using Eq. (89 1 in Eq. (31) and averaging over initial hadron spins gives 






s Xi X 2 



s X x X 2 ij 



d 4 h d A k 2 
(2^)4(2^)4 

d 4 k x d 4 k 2 



|M (l3) | 2 (2^) 4 <5( 4 )(fci + fc2-fc3-fc4) 



(2tt) 4 (2tt) 4 



xU^ul^G^' l G u ' jr Ii 



(k 4 ) H^ CC H[*'i u{kA (2ir) 4 6 (4 Hki + h - k 3 - k 4 ) 



(91) 

The sums over X\ and X 2 represent the sums and integrals over the final states in the C and U subgraphs. The spin 
sums for the incoming hadrons are included in these sums as well. In Eq. (91), the explicit momentum arguments of 



the various factors have been dropped for convenience. The outgoing quark has spin label s. The minus sign from 
the spin sum on the prompt photon has been absorbed into the definition of C. 



To put Eq. (91) into a form closer to Eq. (51), we now apply the minimal partonic approximations of Eqs. (38 
40). Specifically, we neglect the smallest (order ~ A 2 /Q) momenta k± and k^ everywhere except within their own 
respective parent subgraphs, C{P\\k\) and U(P 2) k 2 ). Then the right side of Eq. (91 ) becomes 
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(92) 
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In lines 1-3 we have used the integrals over k2t, k^ , and k 2 to evaluate the 5- functions in Eq. (91 ), and the positions of 
the integrals over k± and k 2 have been arranged into factors corresponding to the respective subgraphs for hadron-1 
and hadron-2, with the factors that correspond to effective TMD PDFs separated by parentheses. The use of the 
kinematical approximations in Eqs. ( 38][40 ) is symbolized by the "hats" on 7?lo- After the second equality in Eq. (92 1, 
we have identified effective tree-level TMD PDFs: 






^ 2 (-P 2 ;fc2(z2,k 2t ; 



x 2 



(Ik 



2n 



Q, 



r1h + x / 



-ft 1 l 



"■1 ^It 

ni ■ k 2 



i v •? 



l l K 2 



2/ 



ni • fc 2 



Q° 



(93) 
(94) 



Also, after the last equality of Eq. ( 92 ) we have restored the explicit momentum arguments to make clear the connection 



with the statement of TMD-factorization in Eq. (37). Equations (94) are consistent with the basic structure of TMD 



PDF operator definitions that are typical in a generalized TMD parton model, and they are consistent with criteria 
of Sect. IV A Hence, in Eq. (|92|) we have recovered a result consistent with TMD-factorization under the maximally 



general criteria. In these graphs there is not yet sensitivity to gauge links type effects; for that we will need to consider 
the extra soft gluons of Fig. [3j 



With the tree level cross section now in the form of Eq. ( 37 ) , one may proceed with the classification of quark and 



gluon polarization dependence following the normal steps in Sect. IV B The contribution to the cross section that 
corresponds to a generic quark polarization projection, T q , and gluon polarization, T g , is written as in Eq. (51): 



da^ r ° 



^EE / ^ <^f* ] (ft;te,k*)) *#f J (PMx^u)) x 



(27T) 
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(95) 



For example, the standard azimuthally symmetric contribution is, from Eqs. (41|49): 



(27T) 
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u(fc 4 )J?LO CT (ki{x 1 ,k lt ),k 2 {x2,cit 



7" ^io.a (*i(a;i, kit), ^2(0:2, 



qt -ki t ),fe3,fc 4 



) u{ki 



(96) 



As a final step, one may average over quark and gluon color to reproduce the standard expressions for the quark and 



gluon TMD PDFs with a color averaged hard part, though this is not necessary to satisfy the criteria of Sect. IV 



To summarize, in this section we have verified that the generalized parton model picture of TMD-factorization 
arises in a detailed tree-level, single-gluon treatment, and we have shown that it is consistent with the maximally 
general statement of TMD-factorization from Sect. IV Tallying the powers of Q in the underbraces of Eq. (95 1 verifies 
that it is a leading-power contribution. 

Later we will follow similar steps when we include one extra gluon from hadron-2, but there they will fail to 
reproduce a factorized form, even under the loose criteria of Sect. |IV A| 



VI. SPECTATOR-SPECTATOR INTERACTIONS 



The next step is to examine the contribution from one extra soft gluon in Fig. [3j If the generalized criteria for TMD- 
factorization are to be respected, then the sum over all such graphs must allow contributions from the extra gluon to 
be factored into separate contributions in the upper and lower subgraphs. The single extra soft gluon contributions 
include both spectator attachments as well as attachments to active partons. In this section we will deal separately 
with the spectator-spectator type graphs, arguing that any leading or superleading contributions cancel. Later we 
will deal with the active parton attachments and find that they lead to TMD-factorization breaking. 
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FIG. 6. The amplitude in Eq. (97 1 with a single spectator-spectator interaction I. 




P, - k, 1\ -k,-l 




Pi - h p.- k, + 1 



(a) 



(b) 



FIG. 7. Spectator-spectator interactions that cancel in the inclusive cross section. The active gluon is a G(k2) gluon which 
yields only leading (non superleading) contributions. (As in Fig. pj] we restrict consideration here to spectator attachments for 
I in the upper and lower bubbles.) 



The relevant amplitude for the spectator-spectator case is shown in Fig. [6j and we extend the notation of Sect. |V A| 
by writing it as 



, f spcc- 



spec 



= U(P 2 ;k 2 ,l) 



a/3;A<i/4 



5/4^1 5/4^2 



0/3:^1^2 



(ki,k 2 ,k 3 ,k4,,l,Pi). 



where 



R 



a/3;fj,iV2 



{k u k 2 , k 3 , k 4 , 1, P t ) = u{k 4 ) ff L o(fci, k 2 , k 3 , k 4 ) a ^° C\P U fci, I) 



/3;M2 



(97) 



(98) 



includes both the bottom bubble £' and the hard subgraph i?LO from Eq. (90). (See the labeling in Fig. [6]) The 
partons are nearly on-shell, k\ ~ k\ ~ I 2 < A 2 , but now the four-momentum may be shared between k 2 and I. As 
usual, we are interested in the region k\ : t % A and k 2) t ^ A. In Fig. R3[ we restrict consideration to graphs with 
spectator attachments for / inside the upper and lower bubbles. 

The superleading contributions come from the contraction of R a l 3 ' + + with 14(P 2 ; k 2 , l) a ^' . The region ^ <Q 
is forbidden because it pushes the quark propagator k\, and propagators with momentum P\ — k\ inside £', far 
off shell. Thus, in the k^ <C Q region, there are at least two extra powers of suppression so that the overall 
contribution is suppressed by at least one power of Q relative to the leading power, even in the contraction of 
R a P'< ++ (ki,k 2 , fc 3 , fc 4 , 1, Pi) with U(P 2 ; k 2 , l) a/3 ' . With k^ ~ Q, both the l + and the l~ components of the i-gluon 
are trapped at order A 2 /Q by initial and final state poles in U{P 2 \k 2 ,l) a ^ ^^ and in £{P\, fei, IY 2 '^ . That is, 
|^ + Z~| -c |1t| 2 % A 2 . (This is the Glauber region discussed in Sect. IIIC ) But then the kinematics of the hard 



subprocess becomes identical to the single gluon case of Sect. |VJ with the only difference being the slightly shifted 
momenta P 2 -~ k 2 — I and Pi — k\ + 1 of the inclusive final state remnants. Therefore, we may once again decompose 
the 3 M ' lAtl in Eq. (97 1 into G(k 2 ) and K{k 2 ) gluons as in Eqs. (76177) to find that the K{k 2 ) terms cancel to leading 
power in an argument identical to that of the previous section. Since k% ~ Q, the remaining G(k 2 ) term leaves only 
leading, not superleading, contributions. 

At this stage, the cancellation of the remaining spectator-spectator Glauber interactions in the inclusive sum over 
the final states of C and hi in the inclusive cross section follows steps familiar from other hadron-hadron processes like 
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the Drell-Yan process. If a single extra spectator-spectator gluon is exchanged, the graphs that contribute at leading 
power to the cross section are as shown in Fig. [71 where the k 2 gluon is labeled by G{k 2 ) to emphasize that these 
graphs are leading and not superleading, according to the argument above, and involve only G-gluon attachments in 
the hard subgraphs. The extra I momentum is routed entirely through the hadron-1 and hadron-2 bubbles, which 
are summed inclusively in the cross section. The cancellation is between the two final state cuts shown in Fig. [7] with 
different final state momenta inside hi and £ . 



VII. ANALYSIS OF TWO GLUONS 



The objective for the next two sections is to show, for the process of Sect. II A| that there is a leading-power 



violation of the loose conditions for TMD-factorization from Sect. IIVI The demonstration will be for the case of a 
single extra soft gluon radiated from hadron-2 and attaching to active partons, with the relevant Feynman graphs at 
the amplitude level shown in Fig. [Hj 

Graphs where one gluon couples to the other before the hard collision were already included in the single gluon 
case of Sect.[Vj The possibility of TMD-factorization violation coming from spectator-spectator type interactions was 
addressed in Scct. |VI| Spectator attachments in hadron-1 attaching to active partons are sensitive to non-perturbative 
structure inside CJ . Therefore, we now focus only on gluons attaching the hadron-2 subgraph to active partons, as 
shown in the graphs in Fig. [51 There are only gluon-fermion attachments in the hard part so ghosts do not contribute 
at this order. 

The steps in this section are nearly the same as those of Ref. [ST] , but are tailored to exhibit the role of intrinsic 
transverse momentum. 



A. Power Counting 



For the two-gluon case, the separation into blocks according to Eq. (32) becomes 



M (2g) = U(P 2 ; l u l 2 ) a ^> 2 V^! 9^ 2 u{k A ) H L {k!, k 2 , k 3) k 4 ; l u l 2 ) a ^^L(P 1 ; k 4 ) 

= u{p 2 - h,i 2 ) afiM v^ 9^ 2 R aP '^ 2 (fci. h, h,k A ; h,k, Pi) ■ (99) 



On the second line we have defined an R a P<viw for the two W-gluon case, analogous to Eq. (69) for the one U- 
gluon case. With an extra soft/collincar gluon now entering the Hi,{kx, k%, k 4 ; l 4 , / 2 ) a/3 ' /ilA ' 2 subgraph, there is another 
fermion propagator in the unfactorized hard subprocess, so each Lorentz component of Hi,{k\, k 3 , k 4 ; l±, l^P'^ 1 ^ 2 now 
acquires an extra power of \/Q relative to the one-gluon case in Eq. ([68]) . The lower quark subgraph, C{Pi;k\), is 



exactly the same as in Sect. W\ so all components of R a P^ 1 > 12 are now of size ~ 1/Q (as usual, up to factors of order 



A and logarithms from renormalization) . By contrast, the one-gluon case in Eq. (71 1 was order ~ Q°. 

The two-gluon subgraph hl a P''^ ltJ " 2 now has two Lorentz indices, and in frame-1 its highest power-law behavior is 
from the hi ~ Q 2 components. The powers for the individual Lorentz components are 

(100) 

(101) 

U*S ~ U+- ~ U— - Q° , (102) 

(103) 
(104) 





U— r 


-Q\ 


u- } 


- W- r 


-Q, 


U + - r- 


-W- + r 


-Q°, 


w+ 


- U +J r 


-i/Q, 




U + + r 


-i/O 2 



The components in Eq. ( 100 ) contribute superleading ~ Q terms graph-by-graph in the amplitude; the U a P'' ~ Q 2 
components of the Zi a/3;Atl/J 2 subgraph multiply R ^'^ 1 ^ 2 ~ \jQ components in Eq. (99). Recall that, while overall 



superleading contributions must cancel in the final analysis, superleading subdiagrams can contribute uncanceled 
leading terms to the amplitude when they multiply subleading factors in the rest of the graph. 
Next we implement the Grammer-Yennie decomposition on Eq. ( |99[ ), writing 

ffMi/ii = ^CiWi + G (h)^ , (105) 

9Av* = K{1 2 )^ 2 + G(l 2 )^ 2 , (106) 
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(a) 



(b) 





(c) 



(d) 





(e) 



(f) 



FIG. 8. Graphs with an extra gluon radiated from hadron-2. 



with the K's and G's defined as in Sect. |IIIB| but now with the four-momenta that correspond to each gluon: 

rtl \ - n "i,/4*i.mi 

< - J A'l,//V 1 /ii — i/^iMl 



<- r i'2J/i^/i 2 — 9li' 2 fJ.2 



m ■ h 

ni,n'j2,iM 2 

m ■ i 



■ '2 



tf(iiW> = 



K(l 



2)fJ,' 2 ll2 



m ■ h 
m -i 2 



(107) 
(108) 
(109) 
(110) 



2G 
Using Eqs. (105 1061 in Eq. (99) separates M( 2g ) into the following set of G/K terms: 



M (2g) =U{P 2] hMT P ^ {G{h) ll , ilXl G{h) ll>2 +G{h) ll , ilXl K{l 2 ) lx ^ + 



K{hUi>MbUi» + ^iW^&WJ R^^z ~ Q ■ ( m ) 

The power-laws indicated with each underbrace denote maximum powers term-by-term. Thus, the largest contribution 
to the amplitude from any individual graph is order ~ Q, which is superleading. As in Sect. [Vl we must verify the 
cancelation of such superleading terms in the sum of all graphs by analyzing each of the G/K combinations. 

B. Basic Expression 

We begin by directly writing down the expression for R"^ 1 ^ from the sum of diagrams in Fig. Q. 

R - ieq9s u( ki , S 4 ) | _ h? _ fc2)2 _ + + Graph a. 



q ' "J lv<* z ' q 

7" 1 (h -t x + m q )r 2 (h -h + m q )rt a tP 

[(fc 4 - h) 2 - m 2 + t0] [(fc 4 - fc 2 ) 2 - m 2 q + *0] 

-^)7' T (^i + /2 + ^g)7 At2 
n 2 +i0] [(fci + Z 2 ) 2 -m 2 

7 CT (# 4 + fe + ™<?)7 M2 (#i + /i + m q )Y ll tH a 



Graph b. 



, y/-4 ri • ■■-gi i v/-i i rz . ■■-ij/ ; Graoh c 

[(fc 4 - ?i) 2 - m 2 + i0] [(fci + h? - m\ + »0] 



Graph d. 
Graph e. 
Graph f. 



[(fc 4 + fc 3 ) 2 - m 2 + iO] [(fci + ?i) 2 - m 2 + iO] 
[(fc 4 + fc 3 ) 2 - m 2 + iO] [Oi + l 2 ) 2 - m 2 + iO] 

i^{k i -t2 + ™ q )r{h + ti + in q ) 1 mh a \ 

[{k 4 ~ l 2 f -m\ + iO] [(fci + ^i) 2 - to 2 + iO] J X 

x£(Pi;fci). (112) 

For ease of reference, we have labeled each term in braces with its corresponding diagram from Fig. [HI Note that we 
have restored the Si argument of the outgoing quark wavefunction. 

C. Glauber regions 

The Glauber regions correspond to intermediate states i n the hard subgraph going on shell, so l 2 and /j - must be 



deformed far into the complex plane, as reviewed in Sect. IIIC to ensure that \l 2 and \l x may treated as order 

~ Q. 

In the two-gluon case of Fig. [HJ only one gluon at a time can be Glauber. If both gluons are Glauber, then 
k 2 << k 2 t, which is forbidden by the kinematical restriction k 2 ~ Q. 

As an example, assume the case where Zj~ is outside the Glauber and consider the integration on 1% ■ For the 
final state ^-interactions in Figs. pfa,f), the deformations away from the ^-Glauber poles are downward, into the 
negative half of the complex plane, while for the initial state ^-interactions in Figs. Mc,e), the deformations are 
upward. Therefore, in Figs. |8 a,f) we apply the downward deformation on l 2 until \l 2 \ ~ Q, and we may replace l 2 



by {l 2 -iO) in Eqs. (108 110) 



[G(l 2 )^ 2 ] F s - g^ r 2 _ iQ 2 > ( 113 ) 

[K{h)^ s =^f^. (114) 

The subscripts "F.S." are to symbolize that these are the G(l 2 ) and K(l 2 ) expressions we will use in graphs with final 
state interactions and downwardly deformed l 2 contours. 
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For the initial state ^-interactions in Figs.p[c,e), we implement the upward deformations on l^ until \l\~\ ~ Q. In 
that case we use 



In +i0 ' 



WW=]i.s. 



n l ,(J,'J 2 ^l 2 
In +10 



(115) 

(116) 



The subscripts "I.S." are to symbolize that these are the G{1 2 ) and K(l 2 ) expressions we will use in graphs with initial 
state interactions and upwardly deformed l 2 contours. 

The requirement is that the signs on the iO's near the l 2 = poles must be such that they do not obstru ct the 
required deformations on l^ away from the Glauber region. Expressions exactly analogous to Eqs. (113 116 1 apply 
for deformations of l\ out of the Glauber region. 

Note that the directions of the deformations are different depending on whether the extra gluon attachments are 
in the initial or final state. So arguments using Glauber deformations need to be applied separately to each graph. 



D. G(h)G{l 2 ) Terms 



Recall from the one-gluon case in Eq. ( 78 ) that the G M ^ factors projected only components of U that were power 



suppressed relative to the largest components. In this subsection, we show that this generalizes to the two-gluon case, 
so the G(l\)G(l 2 ) terms are doubly suppressed relative to the largest terms and contribute only subleading terms to 

M (2 9 )- 



If we apply Eqs. (113 1151, with the deformed contours, to Eq. ( |112 | and recall Eqs. (100-102), then we may 
identify the leading power-law for each contribution from the G(li)G(J,2) projection: 



W(Pi; h, l 2 r^'> GQlJ^MkUi* R (ki,k 2 , h, h; hM,Pxr^ u 



I.S./F.S 



•UiPuhMT™ R{kiMMM;hM,Pi) a 



PM 



0° 



I/O 



■ U{P x ;hMT^ j (1{/QI ± iO))R(h,k 2 , h, fc 4 ; h,l 2 ,Pi) apvl3 

v ' v v ' 

~ ~ I/O 2 



>U{Pi;hMY^- (!i/(£ ±iO))R(k 1 ,k 2 ,k 3 ,k i ;l 1 ,l 2 , Pi) a ^ 3 







I/O 2 



■ U{P x ;hMY^— (l\/(li ± M)M/Q 2 ± iO))R(k!,k 2 , k 3 , h; h,l 2 ,Pi) a0;ij 



o 2 



I/O 3 



1/Q. 



(117) 



Here we have written the largest powers, of order ~ 1/Q, in the contraction with the G(li)G(l 2 ) term. The "I.S./F.S" 
subscript is included to symbolize that we apply the Z-contour deformations upward or downward appropriately in 
each term in accordance with whether the gluon attaches to the initial or final quark. 



The leading powers in Mr 2g ) are of order ~ Q° while the largest powers in Eq. (117) are ~ 1/Q. So the G{l\)G{l 2 ) 
terms are subleading. Notice that the validity of the argument relies on the use of deformed l^ and l 2 contours. 



2.x 



E. K(h) Terms 



Next we write the general expression for the case where one gluon is a A"-gluon. Applying Eqs. ( 114|116 l to 
Eq. (|112|) gives 



[K{l 2 )^ l2 R«^-] 



/F.S 



ie q g 2 s u(/c 4 ,S , 4) 

-,1/2 



h{h ~ U + m q )^ (h -h+ m^tf 3 ^ 



(fc 4 - l 2 ) 2 - rn\ + iO] [(fc 4 - fc 2 ) 2 - m\ + iO] ' Jj - iO 



I/O 2 



7" 1 (# 4 - t x + m,)/ 2 (# 4 -$ 2 + m q )^t a t' 3 1 



(fc 4 - /i) 2 - m 2 + iO] [(fc 4 - M 2 - m- 2 + «0] Z 2 



1/Q 2 



7 Ml (#4 " /l + ™ 9 )7 CT (^i + U + m q )t 2 t a tP 1 



(fc 4 - h) 2 - m 2 + iO] [(fci + Z 2 ) 2 - m 2 + iO] ZJ + JO 



1/Q 2 



rih + h+ mMh + A + m q )^HH a i 

\k 4 + fc 3 ) 2 - m 2 + iO] [(fci + Zi) 2 - m 2 + iO] X Z 2 



1/Q 2 



(a.) 



(b.) 



(c) 



(d.) 



7"(^ 4 + ^3 + ™?)7" 1 (& + / 2 + m q )t 2 t a t^ 



(fc 4 + fc 3 ) 2 - m 2 + iO] [(fci + Z 2 ) 2 - m 2 + iO] Z 2 + iO 

., ' 

~i/Q 2 



(e.) 



t 2 (h - U + ™g)7 g (#i + /l + m.h^t" 1 

"T ~f~^ : , „ ~ .-i f ,. . . „ « rr X 



[(fc 4 - Z 2 ) 2 - m 2 + JO] [(fci + h) 2 - m\ + JO] ij - iO 



1/Q 2 



(£) 



y.n 1 ^C{P 1 -k 1 ) . (118) 

)l/2 



In Figs.p[b,d), l 2 attaches to an internal hard propagator so no intermediate line in the hard subgraph goes on shell 
when l 2 approaches the Glauber region; there are no Glauber poles and so no need for replacements like Eqs. ( 113|116 |. 



F. G{h)K{h) Terms 



The "^i" index in Eq. (118) is ultimately to be contracted with a K(Ji)^> ftl and a G(/i)u'ui m t ne amplitude in 
Eq. (111). We first consider the G(Zi)u' w contraction: 



M, 



( 2 9)lG(h)K(l 2 ) 



ll(P2-JiM a ^' 1 ^G(h)^ 1 K(h)^ >2 R a ^' 



I.S./F.S. 



= ie q g 2 s U(P 2 ;h,l 2 ) a ^ ■" u(fc 4 , S 4 ) X 

Q 1/2 



5/4 Ml 



ni ,^li l|tl \ f 2 (.h - / a + m 9 ) 7 ^ (# 4 - # 2 + mJ'ftPr 



l\ J [{ki-h) 2 -m 2 q + i{)] [(ki~k 2 ) 2 - m\ +i6\ l 2 - iO 



i/Q 2 



VlPl '" '- - »' n ' > r (/c 4 - ?i) 2 - m 2 + iO] [(ft 4 - fe) 2 - to 2 + zO] X ^ 



zr - iO 



3/ii/Ji 



3m i/ 



3m i/ 



i/Q 2 



"i,/i'^i,Mi \ .. 7 Ml (& - /l + "igh'^i + / 2 + m q )t 2 t a t? 
17 -iO 



(fc 4 - ?i) 2 - ™, + *0] [(fci + h) 2 -m 2 q + »0] JJ + iO 



i/Q 2 



"l X^i ,/ii \ 7" (h + h + m q )( 2 {h +k + m q )^nH a 



17+ii) ) [(h + k 3 ) 2 - m q + iO] [(h + h) 2 - m 2 . + iO] I 
ni.M'Ji,^ 7 a (h + h + m q h^{h+t 2 +m q )i 2 t a t^ 



l- 



(k 4 + k 3 ) 2 - m 2 + iO (fci+Z 2 ) 



ra~ + i 



0] lz + iO 



1/g 2 
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(a.) 



(b.) 



(c.) 



(d.) 



(e.) 



ni ,Mi^i ,/ii \ / 2 (^ 4 - /a + ™g)7 CT (^i + /1 + m q )^HH a 
3m>i ~ ! x ^ • ~ - ■--- - • ~ " ~ x 



17 + iO J " [(hi - l 2 ) 2 -m 2 q + iO] [(fei + h) 2 - m 2 + iO] l^ - iO 

^VQ 2 



} 



x£(Pi;fci) 

y v 

~Q 1/2 



(£) 



(119) 



Here we have completed the contraction with U{P 2 \ l\, l 2 ) a ^'^ lfJ - 2 in Eq. (Ill) to get the full G(li)K(l 2 ) contribution 
to the amplitude. 

The next step is to exploit the cascade of cancellations that results when / 2 is substituted with the following 
Feynman replacement identities: 



I2 = -(h - /a - m q) + (#4 ~ m q) 

h = Hi + #1 - m q) - (h - m q) 



'2 — \ki k 2 



- li -m q ) 



tl = (h + h~ m q) - (#1 + /l - m q) 



In graphs (a.) & (f.) 

In graphs (c.) & (e.) 

In graph (b.) 

In graph (d.) 



(120) 
(121) 
(122) 
(123) 



Furthermore, Eq. (119) may be separated into real and imaginary contributions by using the familiar substitution by 

(124) 



distributions convenient for evaluating integrals: 

J^— ( --►P.V.-^-t7rf0 1 - a ), 
"1,2 ' iU *1,2 



P.V.^+f7T^ 2 ), 



»0 



(125) 

'1,2 tu '1,2 

where "P.V." is the symbol for the principal value distribution. Thus, each l/(l7 2 ±i0) eikonal factor is rewritten as a 

sum or difference of two distributions, with the relative signs depending on the direction of the contour deformations. 

For this paper, the main goal is to demonstrate the appearance of a TMD-factorization breaking double spin 

asymmetry with S4 and A3 in the final state. In the squared the amplitude, this arises at lowest order from the cross 
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terms between the graphs of Fig. [8] and the lowest order amplitudes in Fig. [5j With just one extra gluon, it is the 
real terms in Eqs. ( |124|125 ) which yield real terms in the double spin dependent cross section with a violation of the 
minimal TMD-factorization criteria. Therefore, to get the double spin asymmetry we focus attention on the principal 
value contributions^] 

Using Eqs. ^20p23| in Eq. fL19| ), we obtain 



[M, 



(2 S ). 



G(h)K(l 2 ),P.V. 



V v '^ V hJ k 

~ Ql 2 V ^Q '^i/Q 



Y(h -h + m q h°t K 7 CT (^ 4 + # 3 + m q )jPt" 



(k i -k 2 ) 2 -to 2 , -MO (fc 4 + fc 3 ) 2 - m 2 + iO 



VQ 




0(1/Q). (126) 



Here we have combined terms u sing the non-Abelian commutation relation [t Q ,£P] = if a ^ K t K . In addition, we have 
used again the result fr om E qs. (85 86) that the action of l/z 1 — m q on £(Pi; fci) only gives terms suppressed by two 
powers of Q. Equation ( |126[ ) is now almost in the form of the leading order result in Eq. (89), with the leading order 
hard part, Eq. (90), and C(P\; k\) factored away from each other on the second line. The only difference now from 
Eq. (89) is in the details of the factor associated with the upper bubble in the first line of Eq. (126). 



In the contraction with U(P 2 ;li,l 2 ) a P' tl1 • > only the transverse components of ^ are non-suppressed. The plus 
components cancel between the two terms in G(?i)„'„ 1 , and the minus components are suppressed by an extra power 
of Q. Therefore, we may make the replacements l\a' ~* htu,' arL d <jy „ 
leading-power errors and rewrite it as 



g i in Eq. ( 126 1 without introducing any 



[ M V9)] 



G(h)K(h),P.V. 



= -e q9 lf apK u(k^S^U(P 2 ;h,l 2 r^ - (g^ - m y l lt , P P.V.p) P.V.p 

~o 1/2 ' ^ L 




Hence, the G(li)K(l 2 ) terms are leading power. 



x g L o(fci,fc 2 MMT'' P \ C{P i .\h 1 ) 

- I/O ~ Q 1 ' 2 



G. K(h)G(h) Terms 



0{1/Q). (127) 



The treatment of the K(li)G(l 2 ) terms is an exact mirror of the previous subsection, but with the roles of l\ and 
l 2 switched. Since this also involves an interchange of the a and /3 color indices, there is a sign reversal relative to 
Eq. (126) due to the anti-symmetry of / Q ^ K : 



[ M ^)] K{ll 



(h)G(l 2 ),P.V. 



= e q9 if apK u(* 4 , St) u{p 2 - h,h) a ^- •»> <fc p - n x y l 2t , p p.v.- P.V. 



1 



,1/2 



Q 



1 




x g L o(fci,fc 2 , fc 3 ,fc 4 ) K; ^ £(Pi;feQ +Q(1/Q). (128) 
~ VQ ~ Q 1 / 2 



At this order the imaginary parts in Eqs. I |124|l25| l give single spin asymmetries. However, for the one-extra-gluon example there is no 
imaginary contribution that violates maximally generalized TMD-factorization. 
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H. K{h)K{h) Terms 



Extra caution is needed in treating the K(l\)K(l2) contributions because they project terms that are individually 
superleading. To get the K(li)K(l 2 ) terms, we contract Eq. (118) with K(l\)„i ^ to obtain an explicit expression for 
the last term of Eq. (jlllh : 



[ M ^)] K(ll) K {l2) = [u{P2-MMr^'^K{l 1 )^ 1 K{l 2 )^ 2 R^^ 



I.S./F.S. 



= ie q g 2 s U(P 2 ;h 1 h) a ^-^ u{k 4 , S A ) x 
~ Q 2 ~ Q 1/2 



x < — x 



1 t2(h-t2 + m q )t 1 (h-h + m q h' T tH a 



K ' [(h - h) 2 - m 2 q + iO] [(h - k 2 ) 2 - m 2 q + iO] l~ - i{) 



1/Q 2 



/i(#4 - /i + m g )/ 2 (# 4 - # 2 + m q )"/°t a tP 



li-iO [{k i -l 1 ) 2 -m 2 q +i0\[{k i -k 2 ) 2 -m 2 q +iQ] I 



i/Q 2 



I /i(^4 - /i + ™g)7 g (#i + / 2 + m q )l 2 t a tP I 



Zf - iO [(fc 4 - M 2 - m 2 + iO] [(fci + Z 2 ) 2 - m 2 + iO] ' q + iO 

i i a {^ i + h + m q )t2{h + ti + m q )t 1 tn a i 



ir+JO [(fc 4 + fc 3 ) 2 -m2+i0] [(A; 1 + / 1 ) 2 -m2+i0] I 



i/Q 2 



(a.) 



(b.) 



(c.) 



(d.) 



1 7 ff (# 4 + # 3 + m g )/ 1 (^ 1 + / 2 + m g )/ 2 P^ 
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(129) 



Counting the powers of Q written underneath Eq. (129) shows the maximum powers to be superleading. Next, we 
repeat the Feynman identity substitutions of Eqs. ( 120||123 ), again using the result from Eqs. (85 86) that the action 
of fy-y — m q on C(Pi;ki) yields terms suppressed by two powers of Q. In the K(li)K(l 2 ) case above, it is critical that 
dropped terms have two powers of suppression because of the potential for extra leading powers of Q in the cross 
section coming from the multiplication with a superleading U(P 2 ; l\, l 2 ) a "'~ , ~ . 

The principal value contributions are again obtained using Eqs. ( 124|125 ). Thus, after the cascade of cancelations, 
and neglecting power suppressed terms, the principle value K(li)K(l 2 ) term from Eq. (99) becomes 
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The factors of Q in Eq. (130) still combine to give superleading powers, and the equation does not appear to be 



symmetric in l\ and l 2 in spite of it being the K(li)K(l 2 ) contribution. However, Eq. (130) now has the basic 
structure of a leading order result like Eq. (89 1 . Therefore, it may be further simplified with another iteration of the 



G/K Grammer-Yennie decomposition, but now on g p > 



g p » = G{k 2 ) p " + K(k 2 ) p p . 



(131) 



The K{k 2 ) term then vanishes up to power suppressed corrections by an argument identical to that of the one-gluon 
case in Eqs. (81-87). The G(k 2 ) term is then the only remaining contribution and it is 
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After the second equality, we have written the result of contracting l\ p i with the G^fo^^-factor. To write the last 
equality, we have used that k^ = V[ + l 2 and k 2tP = l\ tP + l 2tP . By co nstruction, G{k 2 Y — so the super leading 
powers of Q have canceled — notice after the last equality of Eq. ( 132 ) that the result vanishes exactly for p — " — " . 



We are left after the last equality with a leading-power contribution that is symmetric in a relabeling of l\ and l 2 . As 
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a final step we may also make the leading-power replacement l^ —>■ l^ t , l 2 — > l 2t : 
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Here we have used the notation from Eq. ( 90 ) for the leading order hard part 



I. Complete Amplitude 



Adding Eqs. (127 128 133) gives the complete leading-power principal value contribution to the amplitude: 
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This is the same result found in Eq. (81) of Ref. [81j. It also forms part of the gauge link factor for the ordinary 
collinear gluon PDF. 

The next step is to isolate any terms that can violate the maximally general TMD-factorization criteria of Sect. |IV| 
For this it will be convenient to rewrite Eq. (134) in terms of l\t p and k 2tfl using l 2tp = k 2tp — ht, P '- 
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We now organize the terms in Eq. ( 135 ) according to their compatibility or non-compatibility with the maximally 
general TMD-factorization criteria established in Sects. IV A and Sect. IV D We separate Eq. ( 135 1 into the following 
terms: 
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(137) 
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(139) 
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The J factors are defined as 
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The gluon polarization projections are 
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(147) 



Vu^.vXhtY 1 ^ 1 = -nf nf r it . (148) 

Substituting Eqs.(mfll40} into Eq. pSo) reproduces Eq. ( [T35] ). Note that Eqs. ( fl45J|l46[ ) only project U(P 2l h,l 2 )~ J ~ 
Q compon ents, givin g an overall Q° power when combined with the single eikonal denominators of Eqs. (141 142 1. 
Equations (147 148) project the U(P 2 ; l\, l 2 ) a ^' ~ Q 2 components, giving an overall power of Q° when combined 



with the double eikonal denominators in Eqs._|143 [144) 



Now comparing Eqs. ( |137f|148| ) with Eqs. ( |62| - |65| from Sect. |IVD| verifies that Eqs. ( |137| |138| |l39| ) are consistent 
with TMD-factorization with the maximally general criteria. However, Eq. ( 140 1 is a candidate TMD-factorization 



violating term in the form of Eq. (65). Specifically, the transverse momentum, In, in the gluon polarization projection 
of Eq. ( |148[ ) differs from the transverse momentum entering the hard subprocess. Therefore, we have used "F.V." 
(super)subscripts in Eqs. (140 144 148) to label them as potential TMD-factorization violating contributions. 



J. Candidate TMD-Factorization Violating Terras in the Amplitude 



In Eqs. (137 140) we have categorized each term in the amplitude according to its factorizability, and we have 



isolated the term in Eq. (140) that, according to Sect. IV D can violate the TMD-factorization criteria of Sect. IV A 



Writing out this contribution explicitly in terms of its basic factors gives 
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In the last equality, the amplitude has been written as in Eq. ( 62 ) . The order ~ Q° factor on the first line becomes 
a t( e ffy but now it is exactly in the form of a candidate TMD-factorization breaking contribution like Eq. d65~l). The 



gluon transverse momentum, l lt , which determines the polarization projection in Eq. (1481, differs from the actual 



total transverse gluon momentum, k2t, that enters into the calculation of the leading-order hard part. And it is not 
possible to shift Vu.f.v.O-u) 1 into another factor without Hi, or C then acquiring dependence on li t . 

For an Abelian theory, f a @ K _>. Q in which case contributions like Eq. ( 149 ) vanish. Thus the appearance of 



structures like Eq. ( 149 ) is a consequence of the interplay of color, spin and transverse momentum degrees of freedom 



in a non- Abelian gauge theory. 



K. Reduction To Factorization in Collinear Case 



From the form of Eq. ( 149 ), it is evident that any violation of the criteria of Sect. IV A is specific to observables that 



are sensitive to intrinsic transverse momentum. The violation of generalized TMD-factorization is due to sensitivity 
to the difference between the two intrinsic transverse momenta, li t and k 2 t, outside of U. To have this, the observable 
must be sensitive to intrinsic transverse momenta to begin with. 



More specifically, the "F.V." transverse polarization projection in Eq. (148) appears at leading order in the cross 



section paired with the tree-level projection, G(k 2 ) tJ ' 3 , from Eqs. (77 



89) 



Vu,F.v.(htTG(k 2 y 



(150) 



where to simplify notation here we have only written the transverse i,j indices that couple to the hard part. 

If we were to restrict consideration to observables that are insensitive to a small q t in Fig.[l] i.e. that are azimuthally 
symmetric, then the sum over the i and j components is diagonal. Then we could follow Sect. Ill D of Ref. |81j , 
averaging over transverse gluon polarizations in the hard part while separately tracing over the transverse components 
in Eq. (150) in the upper subgraph, IA, for hadron-2. That is, the combination in Eq. (150) could be replaced with 



^■pu,F.v.(h t yG(k 2 



(151) 



in which case the sums over transverse polarization components, as well as the integrations over l lt and k 2 t, take 
place independently of the rest of the graph. The result is a contribution to the normal integrated collinear gluon 
PDF identified in Sect. V G of Ref. [81j . Hence, factorization of the gluon distribution is maintained so long as we 
restrict to collinear observables. 

In the next section, however, we are specifically interested in obse rvables that are sensitive to the direction of q t , 
and therefore we must account for the off-diagonal parts of Eq. (150). 



VIII. FINAL STATE SPIN DEPENDENCE IN THE PRESENCE OF ONE EXTRA GLUON 



The aim of this section is to verify, by directly counting powers of Q, that terms like Eq. ( 149 1 result in a breakdown 



of TMD-factorization with an extra TMD-factorization- violating spin correlation in the final state. Recall the strategy 
outlined in Sect. |IVC] for identifying spin correlations that are specific to TMD-factorization breaking effects. There 
it was noted that a leading-power correlation between the spin S4 of the final state quark and the helicity A3 of the 
prompt photon does not fit into the standard classification scheme. Thus, the strategy of this section is to isolate 
terms that give contributions from these specific spin correlations in the cross section. 

Such extra asymmetries can only arise from contributions that violate the most general version of the TMD- 
factorization hypothesis, stated in Sect. IV Therefore, we may focus attention on the only term, isolated in Eq. (149), 
that violates those criteria. 



A. TMD-Factorization Breaking Cross Section 



The lowest order contribution to the cross section from the term in Eq. ( 149 1 is obtained by adding it to the tree 
level contribution in Eq. (89), taking the square modulus, and keeping the order g® terms. The steps are exactly 



:-!() 



analogous to those of Sect. VC and we write the details here. Using Eq. (31 ), we find 
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In the second equality, we have written out the explicit factors. In the the third equality, we have changed variables 
from l^ to /c 2 — l\ and we have evaluated the momentum conserving (^-functions using the I2 integration. In the last 
equality we have applied the minimal TMD kinematic approximations of Eqs. ( 38p0 ), as indicated by the hats placed 
on -Hlo: an( i we have shifted the positions of the l^ , /J - , fc 2 and k^ integrals into separate factors reminiscent of the 



step in Eq. (92 1. Also, we have explicitly written the final state photon polarization vectors t a x t\° from Eq. (29) 



Eq. (152) has almost exactly the same structure as the simple tree level result, especially if the factors in braces in 



the last equality arc identified as special types of effective TMD PDFs. Compare with Eq. (92); the only modification 



is that we have isolated the TMD-factorization violating projection from Eq. (150) that comes from a single extra 
soft/collinear gluon. 

A natural strategy would be to attempt to rewrite this as an effective TMD-factorization formula. However, it 
is clear now from the structure of Eq. ( 152 ) that this is impossible without the effective TMD PDFs violating at 
least one of the the minimal criteria for generalized TMD-factorization enumerated in Sect. |IV A| On one hand, if 
the Vu,F.v.0-it) fil>l2 ' 1 ' projection is included inside the definition of an effective gluon TMD PDF then criterion 2.) is 
violated because the projection of Dirac components on the fermion line then depends on li t . On the other hand, if 
Pw.F.vXlit)^ 1 '' 2 ' 1 is included in a redefinition of the hard part or in the definition of an effective quark TMD PDF 
then criteria 1.) or 3.) are violated respectively. So there is no possible rearrangement of factors that allows for a 
factorization formula with the properties of Sect. |IV A| Hence, the maximally general form of TMD-factorization is 
broken. 
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B. Effective Factorization Violating TMD Functions 

In spite of the problems noted above with TMD-factorization, let us push forward with an identification of effective 
TMD PDFs, but now allowing them to violate the maximally general criteria for TMD-factorization from Sect. |IV A] 
Once we have organized factors into a set of effective (but TMD-factorization violating) TMD PDFs, we may then 
directly examine the consequences for power counting and the classification of final state spin dependence. 



We use the explicit expression for ff L ' 0a from Eq. (90) and directly apply the expression for the "F.V." projection 



in Eq. (148). Then, dropping power suppressed terms, the factor inside the Dirac trace on the last line of Eq. (152) 
becomes 



TVf.v. (lit) 



IX-yH^l 



St. 



- (u(k 4 ,S 4 ) fh, £(Pi;fc 1 )£(P 1 ;fc 1 ) 7 ^ 7 ^ u(k 4 ,S 4 ) 



+ u{k 4l S 4 ) la l £(P 1 ;/ci)£(Pi;fc 1 ) 7(r ^7 J ' u(fc 4 ,&0) e^ejf 't K t K ' 
P "-"' v - (1l ' '''" ~ l (u(k 4 ,S 4 ) 7 ^7 l £(Pi;fc 1 )£(Pi;fci)7 J >7 CT ' u(k 4 ,S 4 )) el 3 e* x l't K t K ' 
(u(k 4 , St) fh* £(Pn fei)£(Pi; k^Al 3 u(k 4 , S 4 )) c%^'t K t K ' 



s* 



'1 "1 



St. 



u(k 4 ,S 4 )t lt fij a C(Pi; ki)£{Pv, fa)'? 3 H<r> u(k 4 ,S 4 )- 



F.S.I. 



u(k 4l S 4 ) 7^ f lt C(I\; fei) £(P i; fei) 7CT ^7 j «(**, S 4 ) j e^ejf 't K t K ' 

I.S.I. 

" ! " ' ' u(k 4: S 4 ) la d ( lt C(Pi\ki) £(Pn fci) 7 J >7 CT ' «(**, Si) I el 3 e* x ft K t K ' 



I.S.I. 

J2 (fc 4 , S^tu l^CjP^k^CjP^k^Al 3 u{k 4 , S 4 ) el 3 elft K t K ' . (153) 
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After the second equality, we have included horizontal brackets underneath the equation to indicate where the anoma- 
lous fit transverse momentum contraction is an extra soft initial state attachment (labeled I.S.I.) and where it is an 
extra soft final state attachment (labeled F.S.I.). 

It is in principle possible for soft I.S.I, and F.S.I, attachments to cancel. But any such cancellations cannot be general 
because they depend on the separate non-perturbative details of / lt £(Pi; k\) in the initial state and u(k 4 ,S 4 )t u in 
the final state, including sensitivity to the species of incoming and outgoing hadrons and the quark masses. 

Therefore, without loss of generality, we may focus attention on the initial state interactions. These correspond 
to modified TMD-factorization breaking effective TMD PDFs. An analysis similar to what we give below may be 
applied to the extra final state interactions, but this will give modified effective TMD-factorization breaking TMD 
fragmentation or jet functions. 



It is the second and third terms in Eq. ( 153 1 that contain extra initial state Dirac matrix contractions. These 



correspond to a cross term and an s-channel term respectively in the squared amplitude. We rewrite them by first 



defining the effective TMD functions for the factors in braces in Eq. (152), 
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Note the rf/ 14 in the definition in Eq. (154 1 as compared with Eq. ([93J)p] Equation (154) also has an extra power 
of Q relative to Eq. ( 93 1 and is defined with an overall factor of ljrn q to maintain the same units as Eq. ( 93 1 . 
Equation (155) is just the spin sum for the final state quark, but represented in notation reminiscent of a jet or 



fragmentation function. 



Using Eqs. (153 156) inside Eq. (152) gives 
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In Eq. (157), we have pulled the l/£ outside the integrals using the leading-power approximation t sa (fc 3 — k\) « 
— KiPj feJj~. The subscripts "X — term" and "s — chan" label the cross-term and s-channel contributions. 

The above steps for separating the cross section into factors are very similar to those of Sect. |VC| However, 
Eqs. (157 158) violate the minimal criteria for TMD-factorization from Sect. IV A because of the extra momentum 



arguments of $ / p 1 (-P2; ^2(^2, k2t), lit ) and <E> FJ, lPi;ki(xi,k.it),q,\it)- As explained in the last paragraph of 

Sect. |VIII A| at least one such violation of the minimal criteria is required by the pattern of non-Abelian Ward 

identity cancellations in Sect. |VII| 

For the sake of demonstrating with explicit color factors, let us now take the target hadrons P\ and P 2 to be on-shell 

color triplet quarks, and calculate the contribution for averaged initial and summed final state quark color. Also, let 

(2a) / * \ a P K ';3 ., „ „ , 

us consider the component of $ ,p ( P 2 ; fe2(#2,k2t), lit ) proportional to —if a P R t K t K so that we may define 

an effective colorless gluon TMD function via 
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(159) 



5 An equivalent analysis may be performed with $ and/or l lt organized into different factors. We find the above arrangement to be the 
most convenient for illustrating the basic result. 
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Averaging over initial state triplet c olor and summing over final state triplet color then gives a specific overall color 
factor of (TV 2 - l)/4iV 2 . Equations (|157||158|> become: 
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(161) 



In these equations we have used that CC is diagonal in triplet color; the Tr^ {} denotes a trace over color triplet 
indices. 

To obtain an analogue of Eq. (95 ), we first decompose $ , 3 1 ( P 2 ; ^2(^2, k 2 t), lit ) into its possible transverse index 
structures. The only dependence on transverse indices is from lit and k 2 t so the analogue of Eq. (41) is 



^ 2 (P 2 ;fc 2 (z 2 ,k 2t ),l 






*S fl A Pwl (^;fc 2 (x 2 ,k 2 t),lit) - Q° . (162) 
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There may in general be a separate contribution from each effective gluon TMD PDF labeled with the [fc 2 t] and [lit] 
superscripts. 



We will focus first on the contribution from the [fc 2 t]-term in Eq. ( 162 1. Using an arbitrary Dirac projection T q and 
following by analogy with Eq. (51), Eqs. (160 161) then become 
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FIG. 9. Graphs representing the organization of factors in a.) Eq. (1631 and b.) Eq. (1641. Hermitian conjugate graphs should 
also be included. The "(x)" vertex denotes the anomalous factor of rf/ lt in Eq. (154 1, which cannot be associated with a gauge 
link operator in a TMD PDF for either hadron alone. 
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in accordance with Eq. 

In Eqs. (163 164 1 we have reached the analogues of Eq. (151]), but now in a form that accounts for the extra 



initial state interactions. The effective TMD PDFs now depend on extra momenta external to their parent hadrons 
and therefore violate the minimal requirements for TMD-factorization. The organization of factors in Eqs. ( 153 ) 
through (165) is represented graphically in Fig [9] The special "®" vertex denotes the TMD-factorization anomaly 
factor of dJit in Eq. (154) that enters into the trace of Dirac matrices along the lower active quark line. The extra 



fit ¥" $2t from hadron-2 is trapped by the d on the hadron-1 side of the graph. 



C. Power Counting and TMD-Factorization Violating Spin Asymmetries 



In Eqs. ( 163|164 |, the initial state interaction contribution to the cross section from Eq. (1401 has been written 
in a form analogous to Eq. (51), but with the effective TMD functions now explicitly violating the minimal TMD- 
factorization criteria of Sect. |IV~A) We can now examine whether the exclusion of certain leading-power correlations 
on the basis of parity invariance, as discussed in Sect. |IVC| still applies to the effective TMD PDFs in the TMD- 
factorization breaking scenario of Eqs. (163|164). 
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In fact, a close inspection of the more general multi-gluon effective TMD objects like Eq. (165) hints at a large 



variety of extra angular and spin dependencies. It is beyond the scope of this paper to classify all of them in detail, 
so we will instead illustrate the general existence of such effects by focusing on the specific case, already discussed in 
Sect. [IVC| of a final state correlation between photon helicity A3 and transverse quark spin S4. 



1. Initial State Effective TMD PDF 



Applying the projection in Eq. (59 1 to Eqs. (163 165) gives the effective (but TMD-factorization violating) TMD 
PDF: 



(166) 

Now the effective TMD PDF is a function of q and li t as well as Pi and k\ , in violation of the minimal criteria for 
TMD-factorization from Sect. IV A So Eq. (166[) in general gives a non- vanishing leading-power trace, 



$ (2 9 H76] (p i; fc l(a ; l7 k lt ),<Z,iK) ex ^Tr^/i^A] 
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(167) 



The approximation sign is used on the second line because we have neglected subleading terms in the evaluation of 
the trace. 



2. Effective Hard Subprocess 



163) and Eq. (164). It is the trace for the bottom quark line in 



54. Since our main interest is in the dependence on the final state 



There is one more Dirac trace in each of Eq 
Fig. [9] involving the final state polarizations A3 anc 
transverse polarization for the quark, we keep only the leading power transverse spin contribution from Eq. (155): 

A fc4 / 9 (k 4 , S 4 ) -> 75^4^4/^. 



(168) 



Applying Eq. (60), the Dirac trace for the effective hard subprocess in Eq. (163) for the cross-term contribution 
becomes 



Tr 



4^4f A3 ?A 3 



(169) 



(Recall the definition of h from Eq. ffib.) To get a final state double spin dependence, we must extract the imaginary 
part of Eq. (169) to multiply the i in Eq. (167). For this it is useful to note that for arbitrary four- vectors Vi^ and 



Ui M «2i/Im(e^ 3 e^) = A3 V\ M i> 2 i/Im (e^e* ") . 



Then, the imaginary terms in Eq. ( 169 1 are 
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This trace may be straightforwardly evaluated and the leading-power terms kept. It is helpful to note that 
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where the f^ subscript is a reminder that here the transverse indices i, j are defined with respect to frame-3. Evaluating 
the trace in Eq. ( 171 ) then gives 
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Here we have used the definition of h in Eq. J6). On first line, note that the transverse components of k\ are taken 
with respect to frame-3 and, recalling Eq. (26), the frame-3 transverse components of the incoming hadrons and 
partons are of order 



Q. 



The same steps apply to the s-channel trace of Eq. ( 164 ). The result is obtained by replacing 



ntJ 



"it) 



in Eq. (171 1 . The trace is 
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in Eqs. ( 163 



164) 



Analogous steps also apply if we use the [lu] effective gluon TMD PDF from second term in Eq. ( 162 ) rather than 
the [fc 2 t] one. This contributions differ from the above only in that q t — fc lt gets replaced by l\ t m Eqs. (173 174) and 

(2 S ) [h 



3. Double Final State Spin Asymmetry 



By counting powers of Q in Eqs. 



(162 167 



173 174) inside Eqs. (163 164) we are now able to tally the largest 



powers contributing to a TMD-factorization violating final state double spin asymmetry: 
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(175) 



The imaginary i in Eq. (173) combines with the i in Eq. (171) to give a real contribution to the cross section so the 
graphs in Fig. [9] and their Hermitian conjugates add. 

To summarize, power counting implies a leading contribution to the cross section from a double final state S4-A3 

spin dependence. In Eq. (175), we have explicitly written the A3 (£4) to emphasize the final state spin-helicity 

correlation, and we have explicitly displayed the color factor to emphasize the non-Abelian nature of the effect — the 
spin asymmetry exactly vanishes in the Abelian limit of 7V e — >■ 1 while in QCD it comes with a color factor of 8/9. All 
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173 174) into Eqs. (163 164) have been abbreviated by the 



remaining factors from the substitution of Eqs. ( 162 

O(Q y . This corresponds to a convolution integral of effective non-perturbative multiparton correlation functions. 



Recall that the [75] superscript in Eq. ( 175 ) means that we have applied the pseudo-scalar 75 projection from 
Eq. ( [59] ) to the effective TMD PDF in Eq. ( |154[ ). Without the breakdown in maximally generalized TMD-factorization, 
projections like these would be applied only to ordinary TMD PDFs with the minimal TMD-factorization properties 
from Sect. |IV A| They would therefore be discarded on the grounds that they violate parity. 

To totally characterize the final state angular and spin dependence for Eq. (fTl), all other leading-power projections 
should also be accounted for, including the axial vector and tensor terms from Eq. (44) for the quark subgraph. We 
leave a complete classification of the types of general behavior possible at leading power to future work. We also plan 
to obtain explicit expressions for contributions like Eq. (175), for various initial and final state spin configurations, 



by using specific non-perturbative models for the incoming hadron wave functions. 



IX. SUMMARY AND DISCUSSION 



By directly counting powers of the hard scale for the process of Sect. |II A"} we have shown how extra spin asymmetries 
may arise at leading power due to Ward identity non-cancellations that break TMD-factorization. This is sufficient 
to conclude that TMD-factorization-based classification schemes are too limiting to account for all the qualitatively 
distinct types of non-perturbative correlations that are possible in the limit of fixed energy and large Q > A. The 
(N% — 1)/-/V|? factor in Eq. (175) reflects the special role of non-Abelian gauge invariance in producing such extra spin 
dependence. That the leftover terms arise from an incomplete Ward identity cancellation is indicative of the role of 
TMD-factorization breaking effects in preserving gauge invariance over large distances at leading power in the cross 
section. The essential steps for analyzing Feynman diagrams were already presented in Ref. |81j for the collinear 
case. In this paper, we have argued that certain terms in the Wilson line factors correspond to the breakdown of 
factorization when transitioning from the collinear to TMD cases. 

It should be emphasized that the "extra" soft or collinear gluon attachments here are in the regime of large 
coupling, so we are not actually justified in stopping with a fixed number of gluons. Hence, the failure of the extra 
gluons to factorize in perturbation theory corresponds to the influence of extra non-perturbative effects. The types 
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of interactions are similar to those normally associated with Wilson lines in the definitions of TMD PDFs. However, 
for the process considered in this paper, they cannot be attributed to Wilson line operators for any separate external 
hadron alone. Instead, they describe truly global non-perturbative physical properties of the process. 

The specific example used in Sects. |VH| and |VHI| was chosen to illustrate the existence of TMD-factorization 
breaking correlations in the simplest way possible and with only one extra soft gluon needing to be considered. This 
was sufficient for the current purpose of demonstrating a non-trivial interplay between transverse momentum, spin, 
and color degrees of freedom. However, other extra spin and azimuthal asymmetries are also likely to arise following 



similar mechanisms, especially when going beyond the single extra gluon case. For example, the term in Eq. ( 140 ) 
may contribute to an azimuthal angular dependence in the totally unpolarized cross section. Also, if the imaginary 
parts in Eqs. ( 124 125[) are taken into consideration, then a single final state spin asymmetry may be generated at the 



single extra gluon level. A full classification of spin and angular asymmetries in TMD-factorization breaking scenarios 
is left for future work. 

The interactions discussed here and in earlier work on similar topics (e.g., [T3Hl9t [26 j [90 l l91j) are novel because they 
tend to give rise to effects outside of the common classical intuition that originates in a parton model picture. While 
the general appearance of quantum mechanical correlations is not surprising, these specific effects are distinguishable 
by the fact that they are induced by the kinds of large coupling, long-range interactions that are more commonly 
associated with hadronic binding and confinement in separate hadrons. Even for processes where TMD-factorization 
is valid, physical effects from initial and final state soft and collinear interactions are already predicted in, for example, 
studies of T-odd effects such as the Sivers effect [HI ED- Namely, the Sivers function [52 1 155 ] (a type of TMD PDF) 
is predicted to reverse sign between Drell-Yan and SIDIS due to a reversal in the direction of the gauge link between 
these two processes [26]. (For physical pictures of the sign reversal mechanism see, e.g., Refs. [78] [94].) Interesting 
connections between QCD process dependence and a QCD version of the Aharonov-Bohm effect have been suggested 
in Refs. [23 [SS]. 

The existence of TMD-factorization breaking implies generally that unexpected experimental phenomena are pos- 
sible. In this context, it is worth recalling again the production of transversely polarized A-hyperons in unpolarized 
hadron-hadron collisions (spontaneous A polarization) 52J. Past TMD-based explanations have focused on effects in 
the fragmentation process, such as a polarizing fragmentation function. The production of transversely polarized A- 
hyperons in pp — > A^ + X collisions was described within a hybrid TMD/collinear- factorization formalism in Ref. |97j . 
with the spontaneous transverse polarization being ascribed to the TMD polarizing fragmentation function. It was ar- 
gued in Ref. [98] (assuming universality for TMD fragmentation functions) that this should imply similar A transverse 
polarizations in the e + e~ collisions in the ALEPH experiment at LEP [53], and that the failure to observe such effects 
suggests that spontaneous production of A polarizations in hadron-hadron collisions might also require contributions 
from initial state interactions. (Significant A polarizations also have not been observed in SIDIS measurements 54^-) 

However, the semi-inclusive pp — > A^ + X and e + e~ — > A^ + X processes are not differential in a small qt, and so 
are not strictly appropriate for treatment within a TMD formalism. Moreover, for totally inclusive e + e~ annihilation 
processes, twist-three formalisms [9"9lll00] are so far consistent with measurements of spontaneous A polarization. At 
best, information about the TMDs must be extracted from these measurements indirectly through, for example, hybrid 
formalisms like that of Ref. [97 . Furthermore, the e + e~ measurements have been performed with a hard scale near the 
Z-pole, whereas the measurements of spontaneous transverse A-hyperon polarization in the pp — > A^ + X process has 
a hard scale of only a few GeV. QCD evolution likely leads to a large suppression in the e + e~ measurements relative 
to the hadro-production measurements, and this may be sufficient to explain the discrepancy. In addition, current 
SIDIS measurements may be small due to effects specific to small- ir dynamics. Thus, the status of comparisons of 
spontaneous transverse A polarization between different processes remains inconclusive, and new studies are needed. 
Ideally, investigations of TMD universality should be made by comparing processes like pp — > hadron(jet) + A^ + X 
and e + e~ — > hadron(jet) + A^ + X where a TMD-style description is most appropriate, as explained recently in 



Refs. [1011 I102J . New SIDIS measurement may also be useful, as this is a process where TMD-factorization is 
expected to hold. (The universality of fragmentation functions between e + e~ — ► hadron(jet) + A^ + X and SIDIS 
has been shown in Refs. |103l 1104] .) Future measurements may be possible, and the above general observations 
suggest that spontaneous A polarization studies are ideal for comparing TMD-factorization versus TMD-factorization 
breaking scenarios. See also the discussion in Ref. |105| . Alternative explanations for spontaneous A polarization 
have proposed in Refs. [HE] I106J within the twist-three collinear factorization formalism. It is also interesting to note 
that, in Ref. [107 , a possible inconstancy was found with universal unpolarized TMD A-fragmentation functions in 
pp —> A/A + X measurements as compared with e + e~ — > A/A + X. (See also the discussion in Ref. |108] .) 

Other observables are also likely to be useful for probing the role of TMD-factorization breaking effects. The role of 



three-gluon correlators in a twist-three approach has been recently discussed in Ref. [109J . In Ref. |110j . a proposed 
solution to the Sivers sign discrepancy between Drell-Yan and SIDIS [55] was also given in the twist-three formalism. 
Including a treatment of fragmentation in the description of transverse single spin asymmetries was also found to be 
important within the twist-three formalism in Ref. |lllj . In principle, it should be possible to relate the TMD and 
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higher twist collinear formalisms through integrations over transverse momentum. In this way, it may be possible to 
relate the breakdown of TMD-factorization discussed hear to higher twist collinear correlation functions. Studies of the 
detailed properties of final state jets may also help clarify the dynamics responsible for final state asymmetries. Rather 
than the polarizations of final state hadrons, for instance, one may instead investigate the azimuthal distribution of 
hadrons around a jet axis [112H115] . The "CMS ridge" of Ref. [52] has also been explained in terms of effects that 
are beyond the scope of standard leading twist factorization. These include both Wilson line descriptions [63] and 
color glass condensate descriptions [64 66 . Interestingly, both descriptions incorporate subtle quantum mechanical 
entanglement effects. 

Although a thorough discussion of how TMD-factorization breaking effects might be calculated precisely from first 
principles QCD is beyond the scope of this paper, it is natural to speculate that a form of TMD perturbative QCD, 
with a partonic description of the hard scattering, might be recovered by a rearranged treatment of non-perturbative 
contributions. Rather than having separate TMD PDFs and TMD fragmentation functions, for example, one might 
use a set of non-perturbative functions wherein the constituents of multiple external hadrons remain entangled, but 
a hard partonic factor is separated out. We leave the question of whether such a separation is valid to future 
studies. Already, some versions of soft-collinear-effective-theory (SCET) refrain from identifying operator definitions 



for separate hadrons |116j . Other approaches that are likely to be useful in describing TMD-factorization breaking 
effects are those that are now more commonly associated with the physics of high energy scattering by large nuclei, 
where partons are described as propagating through a large size and/or high density gluon field [117H124"] . 
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